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The  stability  of  small  amplitude  oscillations  in  combustion  chambers  is 
analyzed  for  one—  and  three-dimensional  problems.  In  addition  to  combustion  and 
mass  addition  at  the  boundaries,  residual  combustion  and  the  presence  of  particu¬ 
late  matter  within  the  chamber  are  accounted  for.  The  results  for  the  one¬ 
dimensional  problem  introduce  new  contributions,  to  the  balance  of  acoustic  energy, 
associated  essentially  with  boundary  layer  processes  acting  if  there  is  a  component 
of  acoustical  motion  parallel  to  the  surface.  These  are  incorporated  in  the  general 
three-dimensional  problem,  and  are  shown  to  have  a  significant  influence  on  the  pre¬ 
dicted  stability  of  motions  in  a  rocket  motor. 
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NOMENCLATURE 

Admittance  function  for  pressure  coupling  (Eq.  9.*) 

Speed  of  sound 
Eq.  5.38 

Specific  heat  at  constant  volume 
Eq.  5.35 
Eq.  5.27 
Eq.  5.14 

Internal  energy  per  unit  volume  J  C <AT  of  gases 
Acoustic  energy  density  (Eq.  5.2) 

2 

Stagnation  internal  energy,  e  +  u  /2 

Total  internal  energy  of  particulate  matter  (per  unit  volume  of 
chamber) 

2 

Stagnation  internal  energy  of  particulate  matter,  +  0^*72 
Internal  energy  of  gases  entering  at  the  surface 
Time-averaged  acoustic  energy  (Eq.  7.6) 

Time-averaged  acoustic  energy  (Eq.  7.4) 

Force  of  interaction  between  gases  and  particles  (Eq.  2.8) 

Fq.  4.5 
Eq.  4.16 
Eq.  2.21 

Enthalpy  of  gas  per  unit  volume,  e  +  p/p;  also  Eq.  4.15 
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Stagnation  enthalpy 

Stagnation  enthalpy  of  gases  entering  at  the  surface 
Eq.  2.19 
Eq*  3.8 

Time-averaged  density  of  acoustic  kinetic  energy  (Eq.  7.3) 

Complex  wavenumber,  k  =  (oi-io)/a 

Wavenumber  for  one-dimensional  modes 

Wavenumber  for  three-dimensional  modes 

Length  of  chamber 

Mach  number  of  gases  entering  at  burning  surface,  u^/a 

Vector  mass  flux,  pu,  of  gases 

Mass  flux  of  gas  entering  at  the  burning  surface 

Mass  flux  of  particulate  matter  entering  at  the  burning  surface 

Coordinate  normal  to  the  surface  (Section  6) 

Unit  normal  vector,  positive  outward  from  the  surface 
Eq.  2.17 
Eq.  3.6 
Eq.  4.2 

Time-averaged  potential  energy  (Eq.  7.3) 

Pressure 

Mode  shape  for  one-dimensional  oscillation 

Mode  shape  for  three-dimensional  oscillations 

Heat  release,  per  unit  volume,  in  homogeneous  reactions 

Heat  release,  per  unit  volume  of  chamber,  associated  with 
combustion  of  particulate  matter  (Eq.  2.12) 

Perimeter  of  chamber 


NWC  TP  5349 


R  Gas  constant 

Response  function  for  pressure  coupling  (Eq.  9.2) 
r  Radial  coordinate 

r£  Radius  of  cylindrical  chamber 

Area  of  burning  surface  at  the  end  of  a  chamber 
Cross-sectional  area  of  chamber 
Crcss-sectional  area  of  uniform  chamber 

Area  of  burning  surface  on  the  lateral  surface  of  a  chamber 
Coordinate  parallel  to  the  surface  (Section  6) 

Temperature 
Stagnation  temperature 
Nonisentropic  change  of  temperature 
Time 

Velocity 

Velocity  of  gases  entering  the  chamber,  parallel  to  the  surface 

Velocity  of  gases  normal  to  the  burning  surface 

Velocity  of  particulate  matter  (mass-averaged) 

Velocity  of  particulate  matter  entering  the  chamber  parallel  to 
the  surface 


Combustion  of  particulate  matter  to  gas,  mass  per  unit  volume 
per  second 


6||  Eq.  6-19 
61  Eq.  6.20 


5  „  5  Values  of  6  j|  ,  6j.at  entrance  to  the  nozzle 
nil  n-l 

v. _ _  Roots  of  dJ  (x  r)/dr  =  0  for  r  =  r 

mn  m  mn  c 

p  Density 

Pp  Density  of  particulate  matter,  mass  per  unit  volume  of  chamber 

<f>  Azimuthal  angle 

Gradient  in  direction  normal  to  the  surface 

V ||  Gradient  parallel  to  the  surface 

A  .  Eq.  10. - 

m£  n 

(  )'  Fluctuation 

(  )  Mean  value;  also  time-averaged  value 

(")  Amplitude  of  fluctuation,  without  time  dependence 
(  )  ^  Real  part 

(  Imaginary  part 

(  >  Time  average  (Section  6) 
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1.  INTRODUCTION 


It  has  long  been  recognized  (Ref.  1-4,  for  example)  that  a  proper  description 
of  unstable  motions  in  a  rocket  motor  must  account  for  interactions  between  the 
average  and  fluctuating  flow  fields.  This  is  an  important  feature  not  generally 
included  in  problems  of  classical  acoustics.  The  associated  exchange  of  energy 
may  be  either  a  loss  or  a  gain  for  small  amplitude  motions.  Moreover,  the  con¬ 
tributions  may  be  large  enough  in  some  cases  to  be  dominant  influences  fox  the 
stability  of  oscillations. 

Since  residual  combustion  and  the  presence  of  liquid  and  solid  material  sus¬ 
pended  in  the  gas  phase  are  accounted  for,  the  work  covered  here  is  valid  for  any 
kind  of  combustion  chamber.  However,  most  of  the  discussion  will  be  concerned 
with  unstable  oscillations  in  solid  propellant  rocket  motors;  other  applications 
may  be  extracted  as  special  cases.  For  a  solid  propellant  motor,  the  flow  of  mass 
into  the  chamber  has  some  particularly  interesting  consequences.  This  has  become 
more  apparent  with  recent  work  on  "axial"  or  "longitudinal"  waves  in  a  chamber. 

The  main  reason  is  that  when  there  are  wave  motions  parallel  to  the  surface,  the 
mass  entering  normal  to  the  surface  undergoes  inelastic  processes  of  acceleration 
to  acquire  the  parallel  motions.  As  a  result  of  work  on  the  axial  modes,  earlier 
computations  (Ref.  2  and  2)  of  acoustic  modes  have  been  reexamined;  it  appears 
that  contributions  associated  with  the  mass  flux  at.  the  boundary  were  incompletely 
taken  into  account. 

Since  the  configurations  of  experimental  burners  and  rocket  motors  quite 
naturally  fall  into  several  classes,  the  treatment  here  is  also  split.  At  the 
present  time,  there  is  considerable  interest  in  the  axial  modes,  which  are  best 
handled  in  a  one-dimensional  or  modified  one-dimensional  formulation.  However, 
not  only  are  there  three-dimensional  modes  having  motions  predominantly  parallel 
to  the  boundary,  but  also  the  axial  modes  may  occur  in  very  complicated  configur¬ 
ations.  It  is  therefore  necessary  to  have  at  hand  the  complete  formulation  for 
thr ee-d imens iona 1  prob lems . 

The  emphasis  of  the  present  work  is  mainly  on  linear  problems,  and  particu¬ 
larly  those  aspects  related  to  the  interactions  between  average  and  fluctuating 
quantities.  However,  the  important  influences  of  particles  in  the  flow,  burning 
within  the  volume  of  the  chamber,  and  the  exhaust  nozzle  are  included.  The 
first  part  of  this  report  (Sections  2-4)  cover  the  formulation,  including  lin¬ 
earization  of  three  classes  of  problems:  one-dimensional,  modified  one¬ 
dimensional,  and  the  familiar  three-dimensional  problem.  The  modified  one¬ 
dimensional  problem  takes  into  account,  in  the  unperturbed  problem,  variations 
of  cross-sectional  area  along  the  axis;  the  result  ere  therefore  more  accurate 
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than  those  obtained  from  the  purely  one-dimensional  formulation.  A  restriction 
which  should  be  particularly  noted  is  that  for  all  cases  treated  here,  lineariza¬ 
tion  is  carried  out  according  to  a  limiting  process  discussed  in  Ref.  3;  the 
governing  equations  used  throughout  are  valid  only  if  the  Mach  number  of  the  mean 
flow  is  small.  When  the  Mach  number  becomes  "large"  (i.e.,  greater  than,  roughly, 
0.3),  both  the  mode  structure  and  frequency  are  affected  strongly,  in  respects 
which  cannot  be  ignored. 

For  each  of  the  three  classes  of  problems,  a  wave  equation  is  formed.  Solu¬ 
tion  for  the  case  of  harmonic  motions,  covered  in  Section  5,  produces  formulas  for 
the  complex  wave  number  and  hence  the  growth  constant.  Comparison  of  the  formulas 
for  the  modified  one-dimensional  and  the  three-dimensional  problems  shows  signifi¬ 
cant  differences  which  lead  to  the  principal  new  results  of  the  present  work. 

1.  If  the  interactions,  or  coupling,  between  the  wave  motions  and  surface 
combustion  are  nonisentropic,  then  the  formal  representations  of  the 
coupling  are  different  for  the  two  limiting  cases;  no  acoustic  motion 
parallel  to  the  surface,  and  all  acoustic  motion  parallel  to  the 
surface. 

2.  There  are  energy  losses,  for  the  waves  in  the  chamber,  when  mass 
flows  in  at  an  element  of  surface  where  there  is  a  component  of 
acoustic  motion  parallel  to  the  surface. 

Incorporation  of  these  features  in  the  three-dimensional  problem  is  accomplished 
in  Section  6,  and  leads  to  what  is  here  referred  to  as  the  modified  three- 
dimensional  problem.  The  physical  interpretation  of  these  new  contributions  to 
the  growth  constant — that  is,  to  the  energy  balance  for  acoustic  waves — is 
covered  in  Section  7. 

Perhaps  the  most  important  practical  consequence  of  this  work  is  that  it  is 
clear  now  that  in  general  two  functions  must  be  known  to  characterize  completely 
the  coupling  between  transient  motions  and  combustion.  These  are  conventionally 
introduced  as  the  admittance  and  response  functions;  other  definitions  are  of 
course  possible.  The  physical  origin  of  the  need  for  two  independent  pieces  of 
information  is  that  the  coupling  itself  comprises  two  processes.  One  is  mechani¬ 
cal,  which  may  be  interpreted  grossly  as  p-v  work  associated  with  velocity  or 
mass  fluctuations  of  gas  leaving  the  surface.  The  other  is  due  to  nonisentropic 
temperature  fluctuations  at  the  edge  of  the  combustion  zone.  Since  the  tempera¬ 
ture  fluctuations  in  the  acoustic  field  are  very  closely  isenr.ropic — exactly  so 
within  the  analysis  covered  here — a  nonisentropic  fluctuation  at  the  boundary  will 
generate  a  temperature  or  entropy  wave  which  propagates  through  the  chamber  with 
the  speed  of  a  mean  flow.  The  energy  required  to  generate  such  waves  must  ulti¬ 
mately  be  accounted  for  in  the  energy  balance  of  the  acoustic  waves  within  the 
chamber . 

Several  years  ago,  Cantrell  and  Hart  (Ref.  2)  proposed  that  the  growth  con¬ 
stant  for  acoustic  waves  in  a  chamber,  with  mean  flow  and  combustion  at  the  sur¬ 
face,  be  computed  by  examining  the  total  energy  in  the  chamber.  This  is  clearly 
a  sound  idea,  but  care  is  required  in  the  formulation.  The  analysis  of  Ref.  2  is 
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very  briefly  summarized  in  Section  8  to  show  an  important  restriction:  it  is  not 
possible,  without  modifications  which  are  presently  unknown,  to  determine  the 
effects  of  mass  and  energy  sources  distributed  within  the  volume  of  the  chamber. 
This  seems  to  be  a  consequence  of  treating  the  total  energy  of  the  system.  In 
Section  8.2,  the  linearized  equations  are  used  to  form  an  equation  for  the  balance 
of  acoustic  energy  only.  This  produces  directly  a  formula  for  the  growth  constant 
which  is  identical  with  the  earlier  computation  of  the  complex  wave  number. 

The  last  two  sections  ate  concerned  with  applications  of  the  analysis.  Very 
briefly  in  Section  9  a  simple  configuration  of  the  T-burner  is  examined,  to  indi¬ 
cate  how  experimentally  the  response  and  admittance  functions  might  be  determined. 
Those  functions  must  be  known  for  the  analysis  of  the  stability  of  oscillations  in 
a  motor,  a  subject  treated  in  Section  10.  The  measurements  of  Brownlee  and  Marble 
(Ref.  5),  which  have  been  discussed  in  previous  works,  are  reexamined.  A  conclu¬ 
sion  of  soma  recent  work  by  Perry  (Ref.  7)  was  that  use  of  measurements  taken  in  a 
T-burner,  and  the  analysis  existing  at  that  time-  led  to  significant  differences 
between  the  predicted  and  observed  stability  boundaries.  Although  it  cannot  be  es¬ 
tablished  definitely  at  present,  owing  to  the  lack  of  sufficient  measurements,  it 
appears  that  the  naw  results  obtained  here  may  provide  better  agreement. 

So  far  as  the  formal  analysis  of  oscillations  in  a  combustion  chamber  are  con¬ 
cerned,  the  most  important  subjects  not  covered  here  are  the  influence  of  high 
Mach  number,  of  the  mean  flow,  and  nonlinear  behavior. 
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2.  ONE-DIMENSIONAL  PROBLEM 


This  section  is  based  largely  on  Ref.  4.  The  formulation  is  well-known  in 
fluid  mechanics:  a  one-dimensional  flow  with  mass  addition  at  the  lateral  boun¬ 
dary.  What  is  new  ultimately  concerns  mainly  problems  of  waves  in  such  a  flow, 
with  the  principal  source  of  driving  located  in  the  combustion  zone  at  the  lateral 
boundary. 


2.1.  GOVERNING  EQUATIONS 

The  conservation  equations  are 


mass 

(gas) 

mass 

(particles) 


li{  PSc}  +  4(puSC)  =  J*mbdq  +  WpSc 

w  (pPsc)  +  £  VpV  =  -  wpsc 


(2.1) 

(2.2) 


momentum 


-  r  r  ^  ' 

=  us  Jmbdq+upsJm^P)dq 


(2.3) 


energy 


The  first  two  terms  on  the  right-hand  side  of  Eq.  2.3  represent  momentum  added  to 
the  flow  when  the  axial  speeds  of  the  gases  and  particles  entering  from  the 
lateral  boundary  are  non-zero.  The  two  terms  on  the  right-hand  side  of  Eq.  2.4 
represent  corresponding  contributions  to  the  energy  balance.  Heat  released  by 
chemical  reactions  in  the  gas  phase  is  represented  by  Q;  hence,  e  and  h  represent 
only  thermal  energy  and  enthalpy,  respectively. 
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In  this  section,  variations  of  area  will  be  treated  as  perturbations  from  the 
case  of  a  uniform  chamber,  so  the  equations  for  conservation  of  mass  are  written 


dS 


at  dz 


<Pu)  =  ^-Imbdq  -  sT  ST  + 


9P 


(p), 


p  u  dS 

dq  -  -i-B  ^ 
4  S  ^JaT 
c 


-  w 


(2.5) 

(2.6) 


The  momentum  equation  can  be  rewritten  in  the  form 


■  F-0+u 

(2.7) 

r  9u  3u 

F  —  —  p  — =B  +  p  u  s-B  I 

Up  9t  p  P  9z  -1 

(2.8) 

*  s“  +  (“-up*wp 

(2.9) 

T[,Mu.  S^ill 

O  s 

(2.10) 

Thus,  F  represents  the  force  of  the  particles  on  the  gas,  averaged  over  the  cross- 
section,  and  a  represents  a  momentum  exchange  between  the  chamber  gases  and  two 
sources  of  mass:  the  flow-in  at  the  lateral  boundary,  and  the  mass  of  gas  pro¬ 
duced  by  burning  of  particles.  The  influence  of  differences  between  the  average 
momenta  of  the  channel  flow  and  the  mass  added  at  the  boundary  is  contained  in  u. 

The  energy  equation  can  be  written  for  the  temperature  ir»  the  usual  way, 
first  by  expanding  the  left-hand  side  and  using  the  continuity  in  Eq.  2.1  and 
2.2. 


is2  +  ep-S2  +  -k  <pu)  =  r;  »os-eo,J'mbdq+ c  (ePos-epo)Wl,'dq 


De 


p  Dt 


<P>; 


_(e  -e 
v  o  po 


dS 

)w  +  Q 

;  -r  w  g  dz 

C 
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Now  remove  the  kinetic  energies  from  eQ  and  epQ  by  subtracting  u  times  Eq.  2.3 
and  up  times  Eq.  2.8: 

(h03-e0,K.d,J+  r"  (Vs-epo,=fmbP,dq 


pos  po 

U  V 

dS 

+(up-u)F  +  u(o-n)  +(ei>o-eo)wp-|H_£  +  <Q+Qp> 


C2.ll) 


where  Qp  is  the  energy  released  by  combustion  of  particles: 

QP  =  "  Lpp  “5?  +  Ppup  -g^  ] 


(2.12) 


For  reasons  covered  later,  the  total  energy  and  enthalpy  of  the  gases  issu¬ 
ing  from  the  boundary  will  not  be  assumed  equal  to  the. average  values  in  the 
chamber.  Thus,  while  epos  =  ep,1  hos  ^  hQ  and  eos  ^  eQ.  Define  the  difference 
between  eQs  and  eQ  to  be  Ae0: 


so  that 


e  =  e  +  Ae 
os  o  o 


h  -  e  =  (e  +  Ae  )  +  p/  p  -  e 
os  o  o  o  o 


(2.13) 


=  a  /y  +  AeQ  . 

With  the  equation  of  state,  p  =  pRT,  for  a  perfect  gas,  Eq.  2.11  provides  an 
important  equation  for  the  pressure: 

^  +  yp|H+u||  =  J-  (a^+RATQ )  Jn^dq  +  -A  [(u  U)F+  u(0-U)]  -  -gf  + 

C  v  *  c 


2 

+  #[(epo'^',Wp+<Q+V^  (2'l4) 


3.2.  LINEARIZED  PERTURBATION  EQUATIONS 


For  simplicity,  the  mean  pressure  and  temperature  are  assumed  constant.  All 
quantities  are  written  as  sums  of  mean  values,  (“)»  and  fluctuations,  (  )',  and 
squares  of  fluctuations  are  dropped.2  Then  Eq.  2.7  and  2.14  lead  to 


1  The  kinetic  energy  of  the  particles  is  assumed  to  be  small  compared  with 
the  internal  energy. 

2  Throughout  this  report,  the  primed  (^)*  quantities  contain  dependence  on 
time,  while  later  quantities  signified  by  (  )  do  not.  Thus,  for  linear  harmonic 
motions,  the  fluctuation  of  pressure  is  p’  =  p  exp  (iakt). 
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P'^  +  ’aT  =  "P  g|  <UU')+F'  -a*  +  u« 

3p'  ,  —  3u*  —  3p'  .  du  ,  n. 

-fr +V P-e-  =  -VP'-gj+Pi 

where  Pj  is  the  perturbation  of  the  right-hand  side  of  Eq.  2.14, 

12  o  r  <3S 

Pj  =  j-(a  +RATo)Jmbdq+  ^  [(u  -u)F+u(o-u)l-  ^  -g£  + 
c  v  c 


(2.15) 


(2.16) 


2 

+  '^[(epG-T-,WD+(Q+Q0)]  • 


(2.17) 


The  common  linearized  problems  encounted  in  motors  are  most  conveniently 
analyzed  by  working  with  the  wave  equation  for  pressure  fluctuations.  This  is 
constructed  by  differentiating  Eq.  2.16  with  respect  to  time  and  inserting  Eq. 
2.15  for  3u'/3t  to  give: 

(2.  IS) 

With  ^  * 


hl  =  ~PTZ  (2-19) 

3z  a  a  a 


Solution  to  Eq.  2.18  will  eventually  require  a  boundary  condition  for  p',  set 
according  to  Eq.  2.15, 


with 


=  _f 

3z  1 


fi  =  p  +  p  h  ^u,)  ■  (F,-0,+|j,)  • 


(2.20) 


(2.21) 


In  order  to  obtain  quantitative  results,  the  expressions  for  F',  o',  y*,  and 
P*  are  required.  What  one  uses  for  F'  depends  on  the  law  of  force  between  the 
particles  and  the  gas:  Eq.  2.8  is  merely  the  definition  of  F;  the  force  of  inter¬ 
action  is  required  in  order  to  determine  u„  and  hence  u'  and  F*.  That  problem 
will  not  be  discussed  here.  On  the  other  nand,  o,  and  to  a  lesser  extent  y,  are 
central  to  later  discussion.  It  is  a  straightforward  matter  to  find  o*,  and  y*: 


a'  = 


[u'  Jm^dq  +  u'  JmJP*dq]  +  (u'-u' ) 


i  jw 

P  P 


The  terms 


u'  •  sLC(^)jn7bdq+  (U'  )Jmbdq] 


uJ'mbdq  ■  upImbP)'dq  ■  (u-“p,wp 


(2.22) 


(2.23) 
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and  similar  terms  in  y'  are  dropped  as  higher  order  contributions.  The  reasons 
are  :  and  m^P)’  are  of  order  u,  as  shown  by  computations  of  the  combustion  pro¬ 

cesses  at  the  surface,  and  (u-u  )  is  a  small  quantity  because  the  particles  may  be 
assumed,  usually,  to  be  sufficiently  small  that  in  the  steady  state  they  very  near¬ 
ly  follow  the  gas  motions,  so  that  (u-Up)w'  is  in  practice  a  negligibly  small  term. 


The  terms 


(u  -u)F  »  “(o-yt)  >  (ep0*T^ 


also  produce  fluctuations  of  higher  order  (see  Ref.  4  for  a  discussion  of  the  first 
two)  and  may  safely  be  ignored  in  linear  problems.  In  the  term  a2/y  +  Ae0,  the 
kinetic  energy  produces  a  second  order  term,  and  with  a2  =  yRT,  one  finds  the  per¬ 
turbation 

[(a^+RATQ ) J J 1  =  yR(T+~  ^m^dq+YRtT'  jJ^dq  .  (2.24) 

It  is  assumed  that  the  average  temperature  of  the  gas  leaving  the  surface  is  equal 
to  the  average  chamber  temperature,  so  A?  *=  0,  but  differences  in  the  fluctuations 
will  still  be  allowed;  one  then  has  for  Pj  : 


Pi  =  [T  Jmbdq+(T’  +^y-)j£ibdq]-  -X-  (P>u+pu* 


+r[e;wpv;+(Q+Q^  * 


(2.25) 


The  fluctuation  T'  is  the  isentropic  value  associated  with  acoustic  waves  in  the 
chamber,  while  AT*  is  a  nonisentropic  fluctuation  associated  with  the  unsteady 
combustion  processes  (see  Sections  6  and  8).  This  completes  the  preparations  re¬ 
quired  for  calculations  of  strictly  one-dimensional  problems. 
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3.  MODIFIED  ONE-DIMENSIONAL  PROBLEM 


When  calculations  are  based  on  the  formulation  of  Section  2,  variations  of 
area  are  treated  as  perturbations,  a  point  which  will  become  more  apparent  later 
in  Section  5.1.  It  is  more  accurate  if  the  influence  of  nonuniform  cross-section 
can  be  accounted  for  in  "zeroth  order",  and  in  certain  practically  interesting 
situations  this  can  be  done  without  excessive  difficulty.  To  do  so  requires  only 
relatively  minor  alterations  of  the  formulation  just  covered.  Formally,  the 
difference  consists  mainly  in  retaining  the  area  Sc  on  the  left-hand  side  of  the 
governing  equations. 


3.1.  GOVERNING  EQUATIONS 

The  continuity  equations  for  the  gases  and  particles  are  retained  in  the 
form  in  Eq.  2.1  and  2.2;  the  momentum  in  Eq.  2.2  is  written  (cf.,  Eq.  2.7): 


+  P”S„lr  +  S„|E  =  SIF-o+u)  . 


(3.  1) 


at  *c  dz  '  c  a*  -  c'A  v  * 

Similarly,  the  energy  equation  2.11  is  written  as: 

2 

0  rp  0  m  0  M 

pCvsc  "aF  +  pCvuSc  aF  +P  aF  ("V  =  (T  +  Aeo)Imbdq+sc(V“,F+u(0'u)sc 

+sc»P<V-%>+,Q+Vc  •  (3'2) 

One  then  finds  that,  corresponding  to  Eq.  2.12,  the  equation  for  the  pressure  is 

M  (pScHw  ~(uSc)+uSc  If  =  (a2  +  RaTo)Jmbdq+-i  Sc[(u  -u)F+u(a-U)1 


+  ^Sc«V-T-,WP+(Q+0P,:)  • 

3.2.  LINEARIZED  PERTURBATIOr  EQUATIONS 

The  perturbation  equations  for  velocity  and  pressure  are  notr 

(Scu’)  +  S_-fe-  =  -  pS_  (uu1  HS  (F’-o’+u1 ) 


(3.3) 


c  dz 


PSC  (u u'  HSc(F« -o*+u' ) 


(3.4) 
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-^(p’S^+vp-Jrfu'Sj 


8z 


—  3p’ 

=  -uS 


3z 


•VP'£(>*c>+ScPU 


(3.5) 


where  ScP'n  is  the  fluctuation  of 

Scpu  =  (a2+ATo)J’mbdq+-^-Sc[(u  -ujF+^a-n)]  + 

v  p 


+ 


Hence,  the  wave  equation  for  the  pressure  fluctuations  is  found  to  be 


(3.  6) 


(3.7) 


with 


- 


hu  = 


-  p  2  (uu1 )+  9z  [Sc(F!  -CT'+U  )]  +  _2  ^  s  dz  (aSc)  - 

3z  c  a  a  c 


-  3 


-  P 


3z 


d-f,nS 


3P 


11 


-2 

a 


3t 


(3.8) 


The  associated  boundary  condition  on  p’  is  exactly  (Eq.  2.20)  still,  since  the 
momentum  equation  is  unchanged. 


While  F'  is  still  given  as  the  perturbation  of  the  right-hand  side  of 
Eq.  2.8,  and  o'  is  correctly  given  by  Eq.  2.22,  P’u  differs  from  P{  because 
the  term  ypudfcnSc/dz  in  Eq.  2.15  does  not  appear  in  Eq.  3.6.  Thus,  P'u  is 
given  by 


p*l  =  [  a2 Jm^dq+^RTHRAT'  )J zn^dq^  [(eJwp+epwM+  (QW)  ]  .(3.  9) 


R 
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1 


1 


ft 


|  r  4.  THREE-DIMENSIONAL  PROBLEM:  CLASSICAL  ACOUSTICS 

£  : 


B? 


E 


The  analyses  of  Ref.  2  and  3  which  are,  so  far  as  the  discussion  here  is  con¬ 
cerned,  the  most  general  linear  treatments  available,  are  based  on  modified  forms 
of  classical  acoustics  theory.  The  principal  modifications  are  due  to  the  pres¬ 
ence  of  a  nonuniform  mean  flow.  In  Ref.  3,  the  calculations  were  based  on  inte¬ 
grals,  of  time-averaged  quantities,  over  the  chamber  and  Lts  boundary;  the  de¬ 
velopments  of  Ref.  4  were  founded  on  the  differential  equations.  The  final  re¬ 
sults,  when  comparable,  agree  exactly;  but  for  present  purposes,  it  is  preferable 
to  work  with  the  differential  equations.  Comparisons  of  the  two  approaches  will 
arise  in  subsequent  discussion. 


4.1.  GOVERNING  EQUATIONS 

In  the  differential  equations  for  three-dimensional  problems,  no  boundary 
sources  appear  explicitly,  so  one  has  for  the  conservation  laws,  instead  of  Eq. 
2. 1-2. 4, 


mass 

(gas) 

V-  (pu)  =  wp 

(4.  1) 

mass 

(particles) 

9p 

at  rpP  p 

(4.2) 

momentum 

(pu+ppup) +v-(puu+pupup)+vp  =  0 

(4.3) 

"nergy 

0  *  ■> 

-37  (pe  +p  e  )+V*  (pue  +p  u  e  )  +  V.  (pu)  =  Q 
at  o  rp  po  o  *p  p  po  ’ 

(4.4) 

The  momentum  equation  for  the  gas  phase  is 


where 


p^i  +  pu-  Vu  +  7P  =  F  -  a 


r  ou  -i 

F  =  -  I  p  -3^  +  p  u  -  Vu 

LKp  dt  p  p  p J 


(4.5) 


(4.6) 
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I 


m 


m 


and  a  =  (u  -  u  )w  (4.7) 

P  P 

Consequently,  the  energy  equation  with  temperature  as  the  independent  variable  is 

3  T  - *  — ♦  *4  — ♦  —4  >4 

pc  or  +  pC  u.7T+pV-u  =  (u  -U).  F+u.  0+(e  -e  )w  +(Q+Q  )  (4.8) 

v  oi  v  p  po  op  p 

and  the  equation  for  the  pressure  is 

+  ypV.  u+u- Vp  =  t(Up-u)- F+u*  u+(ep0-eo)Wp+ (Q+Qp)l+RTw^  (4.9) 

Variations  of  cross-sectional  area  of  course  do  not  appear  explicitly;  they  enter 
the  problem  when  the  three-dimensional  differential  equations  are  solved  in  a 
specified  geometry. 


4.2.  LINEARIZED  PERTURBATION  EQUATIONS 

In  Ref.  3  the  limiting  process  for  finding  the  linearized  equations  has  been 
described.  This  will  produce  the  equations  usually  used,  in  which  terms  of  first 
order  in  the  wave  amplitude  and  mean  flow  Mach  number  are  retained.  The  sets  of 
equations  already  found  in  Sections  2  and  3  are  essentially  based  on  that  proce¬ 
dure,  although  they  have  not  been  so  rigorously  deduced  here. 

Prom  Eq.  4.5  and  4.9,  one  finds  the  linearized  equations  for  u'  and  p’: 

—  9~*'  —  Zl  _  _.  _ 

p  +  vp'  =  -p( it*  Vu'  +  u1  •  tfu  )+  F'  -  CT* 

a  i  _  -  rS  - 

~£-  +  "ypV •  u*  =  -u*  Vp'.-yp'V-  u  +  P' 

with 

p.  a  (V  W?  +  SpWp  +  (Q'+Qp)]* 

The  linearized  wave  equation  for  pressure  fluctuations  is 

.  h 

a  dr 


3  Note  that  the  assumption  of  incompressible  mean  flow,  V  »u  =  0,  is  not  made. 
This  is  obviously  necessary  if  there  is  residual  combustion  (wp  y  0) ,  but  it  is 
also  a  formal  convenience  for  later  comparison  with  the  one-dimensional  results. 


(4.  10) 
(4.  ID3 

(4.  12) 

(4.  13) 
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subject  to  the  boundary  condition,  deduced  from  Eq.  4.10,  for  the  normal  gradient 
at  the  boundary: 

n  *  Vp*  =  -f  (4.14) 

The  functions  h  and  f  are  now 

h  =  -pV.  ( u •  7u'+u'-  Vu  )  +  u*  7  +7*  (F1  -0' )+  ^2  OIL  7.  u  _  -i.  ,  (4.  15  ) 

a  a  a 


f  =  5+ p(u.  Vu'+u' •  7u)-  n-  (F*-CT«)- n  .  (4.16) 

As  formulated  here,  the  three-dimensional  problem  does  not  account  totally 
for  all  interactions  between  the  wave  motions  and  the  flow  issuing  from  the  sur¬ 
face.  This  can  be  seen  most  easily  by  examining  the  special  case  of  one-dimensional 
motions.  The  left-hand  side  of  Eq.  4.13  is  then  identical  to  the  left-hand  side 
of  Eq.  2.18  but  the  right-hand  sides  and  the  boundary  conditions  in  Eq.  2.20  and 
4.14  do  not  agree. 


With  variations  only  in  the  z-direction,  Eq.  4.15,  4.16  and  4.12  give 


h  = 


2 


v  9p'  du  _1_  ap- 

-2  at  dz  "  —2  at 

a  a 


f  =  0 


P*  =  [e' w  +e  w4(Q'+Q' )] 

C  p  p  P  p  p 


These  will  obviously  not  reproduce  the  one-dimensional  problem  discussed  in 
Section  2.  Indeed,  there  is  no  possibility  here  of  representing  any  coupling 
with  the  burning  surface!  The  reason,  of  course,  is  that  in  inviscid  classical 
acoustics  the  only  boundary  condition  is  placed  on  the  velocity  fluctuation 
normal  to  the  boundary — this  is  simply  zero,  always,  for  one-dimensional  problems. 
It  is  mainly  this  difficulty  which  motivated  the  work  of  Ref.  4.  Modifications 
of  the  inviscid  three-dimensional  problem  consist  of  approximations  to  the  pro¬ 
cesses  occurring  in  boundary  layers  present,  in  the  real  problem. 

However,  the  case  of  truly  three-dimensional  motions  with  significant  com¬ 
ponents  parallel  to  the  boundary  surface  has  not  yet  been  treated.  This  is  re¬ 
ferred  to  here  as  the  "modified  three-dimensional  problem";  it  is  treated  in 
Section  6.  As  an  aid  to  that  discussion,  it  is  useful  first  to  examine  some  more 
detailed  results  for  the  stability  of  harmonic  motions. 
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5.  CALCULATION  OF  THE  COMPLEX  WAVE  NUMBER 
FOR  THE  LINEARIZED  NATURAL  MODES 


The  normal  or  natural  modes  for  a  linear  system  are  characterized  by  harmonic 
time  dependence;  all  fluctuations  vary  as  exp(iakt)  where 

ak  =  U)  -  ia  (5.1) 

is  the  complex  wave  number.  If  the  growth  constant  a  is  positive,  then  the  ampli¬ 
tude  of  the  waves  grows  with  time.  Now  for  a  linear  acoustic  wave,  the  energy 
density  is  the  sum  of  potential  and  kinetic  energies: 

e  =  +|p(6u)2.  (5.2) 

a  2p  aT 

Here,  <Sp,  6u  are  written  since  they  must  be  real  quantities  to  make  the  energy 
real.4 

All  phases  in  time  are  referred  to  the  pressure  oscillation;  for  harmonic 
motions ,  then 


.  *  at  iiut 

p1  =  pe  e 

(5.3a) 

~x  at  i(t«t+<j>) 
u'  =  ue  e 

(5.3b) 

where  $  is  the  phase  of  the^velocity  fluctuation  relative  to  the  pressure  fluc¬ 
tuation.  In  general,  both  p  and  u  may  be  complex  functions  of  position,  but 
suppose  for  the  present  that  they  are  real.  Then  the  real  parts  of  Eq.  5.3a,  b 
may  be  used  for  6p  and  Su  in  Eq.  5.2  to  give 

e  =  -^2  cos^iut  +  p  (u  )2cos^(U)t+4>)je^0t*'  (5.4) 

a  pa 

and  the  time-averaged  value  is  approximately 


4  One  could  of  course  introduce  complex  conjugates  and  write 

©  =  *  ~  1  ^  p'p1*  +  }  p  u1  •  u»*  . 

2p  a 


14 


NWCTP5349 


ea  =  T'  ^Z-  +  pu2)e2at  .  (5.5) 

pa 

The  time  average  here  is  over  one  cycle  of  the  oscillation  and  is  valid  only  if 
the  amplitude  of  the  wave  doesn't  vary  rapidly:  a  «  w.  It  follows  from  Eq.  5.5 
that  the  growth  constant  is  related  to  the  rate  of  change  of  energy  density  ac¬ 
cording  to 

,  die 

0  la 

2a-  =  —  "dF  •  (5.6) 

e 

a 

For  interpretation  of  later  results,  it  should  be  noted  that  Eq.  5.6  can  be  inte¬ 
grated  over  the  entire  volume  of  the  cavity,  showing  that  2a  equals  the  fractional 
cate  of  change  of  the  total  time-averaged  acoustic  energy  in  the  cavity. 

The  linear  stability  of  the  normal  modes  depends  on  a.  It  is  easy,  following 
the  development  first  given  in  Ref.  3,  to  find  a  formula  for  k2  which  can  be  eval¬ 
uated  using  the  classical  acoustic  mode  shapes.  The  wave  number  for  each  mode  is 
known  to  first  order  in  the  perturbations  (i.e.,  to  first  order  in  the  mean  flow 
Mach  number)  when  the  modes  are  known  to  zeroth  order.  Those  are  the  modes  for  a 
closed  chamber  with  no  combustion  and  having  rigid  walls.  This  is  a  well-known 
characteristic  of  perturbation  expansions  of  the  sort  used  here. 

On  the  other  hand,  the  actual  mode  shapes  in  the  chamber  will  be  distorted 
from  the  classical  modes  by  terms  of  first  order  in  the  mean  flow  Mach  number. 
Moreover,  they  will  be  complex  quantities  of  position,  which  means  that  the  actual 
pressure  and  velocity  fluctuations  will  have  phases  varying  with  position  in  the 
chamber.  All  phases  are  measured  relative  to  the  unperturbed  (classical)  acoustic 
pressure  oscillation.  Thus,  for  example,  Eq.  5.3a  then  becomes 


n.  i(u)t+«j>  ) 
P'=lp|eate  P 


with  <|>p  a  function  of  position.  The  actual  energy  of  the  wave,  Eq.  5.2,  will  also 
depend  to  first  order  on  the  mean  flow  Mach  number.  These  refinements  are  of  no 
consequence  for  the  first  approximation  to  linear  stability. 


5.1.  ONE-DIMENSIONAL  PROBLEM 


With  the  exponential  time  dependence,  all  fluctuations  can  be  written  in  the 
form  p'  =  p  exp(iakt),  for  example.  Hence,  Eq.  2.18  becomes 


(5.7) 


with  the  boundary  condition  at  the  ends  set  by  Eq.  2.20: 
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d£  _  . 

dz  "  ‘  *1  *  (5.8) 


In  the  absence  of  combustion,  Eq.  5.7  and  5.8  become 


i  dS 

.  ,  —  A  1  C 


iE  »  0 

dz 


(z  =  0,L) 


(5.9) 

(5.10) 


Equation  5.9  exhibits  explicitly  the  contribution  of  nonuniform  area  alone, 
when  there  is  no  combustion  or  average  flow. 


The  procedure  for  finding  an  expression  for  the  wave  number  involves  essen¬ 
tially  a  comparison  of  the  problem  of  interest,  including  perturbations,  with 
the  problem  having  no  perturbations.  The  latter  is  simply  the  classical  acoustic 
problem  for  a  straight  tube,  governed  by  the  differential  equations 


Pl  .  j2  *  _o 

7T  l  Pl~ 


(5.11) 


dP, 

=  0  (z=0,L)  (5.12) 


n  A  A 

To  find  the  first  order  formula  for  kz,  multiply  Eq.  5.7  by  p^,  Eq.  5.11  by  p, 
substract  the  equations  and  integrate  over  0  z  <.  L.  After  use  of  the  boundary 
conditions  Eq.  5.8  and  5.12  one  finds 


where 


(5.13) 


(5.14) 


All  terms  in  h|  and  f  i  contain  the  Mach  number  of  the  mean  flow  as  a  factor  (see 
Ref.  3  and  4  for  further  details).  That  is,  all  perturbations  of  the  classical 
problem  are  small,  as  necessary  for  the  calculation  to  be  valid.  Hence,  it  is 
Sufficiently  accurate  (i.e.,  to  first  order)  to  replace  p,  u  everywhere  in 
hj,  fi  and  e|  by  p^  and  u^.  The  right-hand  side  of  Eq.  5.13  may  therefore  be 
evaluated.  To  do  so,  it  will  be  necessary  to  use  the  relation  for  the  acoustic 
velocity. 
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A 


(5.15) 


With  (2.  19)  for  hj  =  hj  exp(iicat)  «  hj  e^i(ik^at)  , 

.i^tp  p  iz  -  f(F  ^a'+u  )-^<^  +  !l^(F  -5  +a  >PJ 

»  0  1  ^  0  CC  0 


(5.16) 


and  with  (2.  21)  for  f^  =  exp(ik^at)  . 


(5.17) 

The  area  Sbe  has  been  introduced  as  the  area  at  the  end  of  the  chamber  over  which 
V  differs  from  zero;  this  may  occur  usually  at  a  burning  surface  or  at  the  en¬ 
trance  to  a  nozzle.  Addition  of  Eq,  5.16  and  5.17  produces,  after  some  further 
integration  by  parts: 


— a<£  _  T 

UP*  Su„  -,L 


z,  2  2  2  . — ,  r  ^  "be  T\ ..\p.2du 

pa  co  f*  a  0 


dz 


L  dp,  k.  L, 

-J  (F  --7  +$')-—  dz"i~  I  PlP*dz 
0  a  0 


(5.18) 


With  Eq.  2.20  and  2-21  one  finds 


L  dp 
“a 

a  -  ~  dz  = 
dz 


i 


L  dP>  2  j_ 

f»  /  •  1* —  .  dz 


pak^O 


dh  r=  (p). 


‘b-  s 


c  0 


.  ,  Jm.  dq  t.- 

p  dz  d  o  ^  S 


L  dp. 

+  f  (u  -u  )  w  dz 
5  p  p  dz 


(5.19) 
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l*Adz  "iak4  PjS,dq|^-i(V-l)Y  l  n  IT-  Jmbdqdz 

'  r* 


L  dp,  dS  , 

1  r  f,  _ c  dz  , 

2  *L  dz  dz  S 

0  c 


k,  L  ,  dS„  , 
t  p  c  dz 

—  J  P^-jrr 

a  0  -  dz  sc 


-i  ~  R  /at  P  Jm^dq 


a 

k 


—  dz 
S 


‘^CT  I  P^[{apwp+  epW^+^'+Q^)]  dz 


(5.20) 


The  -1  part  of  (y-1)  times  the  second  integral  is  combined  with  the  first  inte¬ 
gral  in  Eq.  5.19,  while  the  y  part  is  rewritten  by  use  of  the  continuity  equation 
for  steady  flow: 

_1_  f-  .  du+  u  % 

J  +  s“  "dT  "  "=? 

pSc  c  p 

Then  the  formula  for  k2  has  the  final  form 


© 


(k2-kL)EL  =  {iPakf,[CGP-(i  +  '^')s^  ] 

p  a  co  Q  0  c 


© 


a  p  0 


c  0 


dz  dz  S 


} 


♦{k&=bw*i:-r*2-} 


hi) 

I  *dp< 

c  0 

ivY 


(5.21) 


4  {  V(VVVV+(Q+V]ir  +i~  f^wpdB+f(ia-fip)wp"aF  dzl 


P  'Sc  pa  5*^  P  O’  P  P 
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(5.22) 


The  various  contributions  have  been  grouped  as  shown  according  to  the  following 
general  interpretations: 

©  These  contain  all  terms  representing  the  response  of  burning 
surfaces  (u  is  related  to  the  admittance  function  and  1S5  to 
the  response  function).  In  addition,  the  influence  of  a  nozzle 
at  the  end  of  the  chamber  appears  through  its  admittance  func¬ 
tion  (the  u  term).  The  term  up|  represents  the  convection  of 
acoustic  energy,  by  the  mean  flow,  through  the  end  walls.  It  is 
an  energy  gain  to  the  wave  system  because  it  represents  the  addi¬ 
tion  of  energy  by  the  inward  flow  after  the  flow  has  acquired 
the  acoustical  motions. 


in  These  represent  interactions  between  the  mean  gas  flow  and  the 
acoustic  field.  The  first  term  represents  an  energy  loss  for 
the  waves  in  the  chamber.  It  is  due  to  the  inelastic  process 
by  which  the  incoming  flow  acquires  acoustic  energy.  Note  that 
p|  +  kT2  (dp^/dz)2  is  proportional  to  the  volume  density  of 
acoustic  energy.  The  second  term  is  the  net  (averaged)  work  done 
by  the  acoustic  pressure  against  the  mean  flow.  The  third  term 
shifts  the  natural  frequency.  Its  interpretation  is  not  obvious 
it  arises  from  the  term  yp'u  dfcnSc/dz  in  Eq.  2.25. 


The  first  term  is  associated  with  the  addition  of  particles  from 
the  burning  surface,  and  the  fact  that  they  must  acquire  the 
local  fluctuating  motions.  The  second  term  represents  the  influ¬ 
ence  of  particles  distributed  in  the  volume  of  the  chamber  and 
hence  contains,  partly,  the  attenuation  of  acoustic  waves  by  the 
particle/gas  interactions. 


TVJ.  All  these  terms  vanish  if  there  is  no  residual  burning  within  the 
volume.  Those  containing  w„  or  Wp  as  a  factor  are  present  only 
when  there  is  burning  of  solid  (or  liquid)  particles-  The  first 
terra  is  associated  with  energy  release  and  the  second  term  repre¬ 
sents  a  source  of  acoustic  energy  due  to  the  distributed  source 
of  mass. 


This  represents  the  effect  of  temperature  fluctuations  at  the 
combustion  zone,  which  are  nonisentropic .  Hence,  It  arises  due 
to  the  production  of  entropy,  ai.d  part  of  it  is  an  energy  loss 
for  the  acoustic  waves. 
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The  last  terra  represents  the  effect  of  momentum  fluctuations  at 
the  boundary  layer  when  the  axial  components  of  momenta  of  the 
gases  and  particles  entering  differ  from  those  of  the  gases  and 
particles  in  the  bulk  channel  flow. 


Some  of  the  above  will  be  interpreted  in  more  detail  later. 


5.2.  MODIFIED  ONE-DIMENSIONAL  PROBLEM 

The  calculation  of  the  wave  number  goes  through  in  the  same  way:  the  only 
major  difference  is  that  the  unperturbed  or  zeroth  order  equation  must  now  corre¬ 
spond  to  Eq.  3.7,  which,  for  harmonic  motions,  is 


1 

S 

c 


A. 

dz 


(S 


c  dz 


)  +  k2p 


(5.23) 


The  boundary  condition  is  still  Eq.  5.8,  but  the  unperturbed  problem  is  represented 
by 


c 


(5.24) 


dp 


^  =  ° 

dz 


(5.25) 


Now  multiply  Eq.  5.22  by  Scp^,  Eq.  5.23  by  Scp  and  follow  the  prescription  out¬ 
lined  above  to  find 


,  L 


k  ■  kt+xiIJiAscdsttfAscy 


where  now 


EU  '° 


TT2  -  rL-2o  ^ 

Eu  ~  1  p*  scdz  * 


o 

The  definitions  in  Eq.  2.21  and  3.8  of  f^  and  hn  lead  to 

— a2  L 


(k2-k>U  =  S  Kbel  +iW-l>  ^ 

pa  u  a  u 


L  dp,  k,  L  * 

-  r  (#-s+s)  scdZ-i  r  p  p  s  dz 

0  a  0 


(5.26) 


(5.27) 


(5.28) 
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The  integral  over  <s  is  like  £q.  5.19  except  that  dz/Sc  is  replaced  by  dz; 
that  involving  Pjj  is 


-i  —  J  Pnp.S  dz  =  -iak  J P .Jmhdqdz - *(Y ~l ) Z~  J  pJ  JayM2 
a  0  “  c  b  pa  0  4 


-i  —  fp.  [ e  w  +"e  w  +(Q' +Q)]dz 

V  p  p  p  p  p 


a  Cv  (J 


k  L  * 

-i  — Rj  AT  pt  Jrr^dqdz 
a  0 


(5.29) 


Thus,  Eq.  5.28  becomes 


®  © 

— *2  L  L 

=  {ip^l(fipt^)sbej  -iaV.{,PtJ'Vq,U} 

pa  u  u 

©  @ 

if  tPi+4 J5vjqdz '  i^:«l!)p/r-4<sc",sc'fa} 
a  p  0  k*  ©  a 

+  {  -It  J5ib(p  J  F  Scdz} 


^-‘rfJySpVWn-c 


k.V  L  y  ij  _  't 

5+6)1 S  dz+i-^  J  w  P#S  dz+j  (u-ii  )w  Scdz| 


a  v  0 


p  a  0 


prfc 


P  P 


© 


-i  _i  R  f  ATp^  J nylqdz  -  J  U  Scdz 


(5.30) 


a  0 


The  interpretations  of  the  various  terms  are  the  same  as  those  listed  above. 
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5.3.  THREE-DIMENSIONAL  PROBLEM 


The  calculation  now  centers  on  solution  to  Eq.  4.13.  The  unperturbed  problem 
is  described  by 


_2*  ,  .  2~ 
»  P-.T  + 


■N 


ktpN  *  0 


n.  =  a 


where  N  stands  for  three  indices,  one  for  each  dimension. 


(5.31a) 

(5.31b) 


For  harmonic  motions,  Eq.  4.13  is 


with  the  boundary  condition 


2*  2*  * 

V  p  +  k  p  =  h 


A  -.A  _ 

n  -  Vp  =  -f 


(5.32) 


(5.33) 


Now  multiply  Eq.  5.32  by  p  ,  Eq.  5.31a  by  p,  subtract,  and  integrate  over  the 
volume  of  the  chamber.  After  use  of  Green's  theorem  and  the  boundary  conditions 
in  Eq.  5.31b  and  5.33,  one  finds 

k2  =  k2N  +  {  JpNSdV+  JJ  £  pNdS}  . 


(5.34) 


Again,  of  course,  this  result  is  valid  only  to  first  order  in  the  Mach  number  of 
the  mean  flow.  Thus,  in  fi  and  f  one  must  use  the  zeroth  order  approximation,  p  , 
to  the  mode  shape  (cf.  remarks  following  Eq.  5.6).  Here,  E^  is 

EN  =  JPNdV  *  <5.35) 


The  acoustic  velocity  is  now  given  by  the  formula 


"N 


pakN 


7pN  • 


(5.36) 


Hence,  with  n  given  by  Eq.  4.15,  and  f  by  Eq.  4.16, 

~  J*u.  VpJTdV+2  ^  P-2”  “ 

2a 


j-fipNdv  =-~Ipn^c dV+i J u.  Vp2dv+i 2L  Jp2V.udV+JpNV.  (iJ )dv 

(5.37) 


J  ^PNdv 


£?  =  ip  akN  ^PNdS  +  IT—  CpNdS  -  §  n-  (F-0)pNdS 

akN 
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where 


C  =  u  ‘  V(VpN)  +  Vpjj  •  u 


(5.38) 


All  terms  In  C  can  be  combined  to  give  an  interesting  simpler  result  as  follows; 
since  in  Eq.  5.34  the  integrals  over  h  and  f  must  be  summed,  one  encounters  the 
following  combination: 

-  I  PN*  CdV+  &&■  CpNdS 
a^N  akN 

=  .  -i-  J7.  (pNC)dV+ri-  J C-  VpNdV+  -±-  $  £.  CpNdS 


1  "  “n5‘  idS  +  ^-  J  C-  VpNd V  +  -i-  §  £•  CpNdS  . 


The  first  and  last  terms  cancel,  while  the  second  one  can  be  rewritten* 

C-7PN  =  Cu-  v(7PN)+vPN'v«i]-  VPN 

=  [V(u.  VpN)-7pNxvx^].VpN 

— ♦ 

=  VPN*  *(«•  VPN) 


Consequently, 


*  7-H"-7pn^n:|-(''-7Pn,v2Pn 
■  7-r^-7PN)vpN]+kN2pN^-7pN 


=-*-  ! 5-  7PN-iv  =  J v  ■ •  [(£  vpN)7pN]dV + i  -2L  j£  7pN2dv 


=  JJ(“-7PN)7PN-fldsti^J^7PNdv 

2a 
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The  first  integral  vanishes  because  of  the  boundary  condition  in  Eq.  5.31b,  while 
the  second  is.  exactly  equal  to  the  second  term  in  Eq.  5.37.  Then  since  the  mean 
flow  field  is  not  assumed  to  be  free  of  sources. 


Hence,  after  a  few  other  simple  manipulations,  the  formula  in  Eq.  5.34  for  the  com¬ 
plex  wave  number  can  be  written  explicitly  as 


© 


—  *2 
uPr 


fc2-l4)EN  =  {ip  =3  ]*  5dS  - 


Q 


{— - - i^Jp^dv} 


pa 

© 


-} 


f 

n 


© 


in 

Jf  .  vpNdv 


+  J V  VV  VV +(°+ V]dV+i  =■  -f  PnVV 

av  r  IT  r  p  a  r 


+Jwp(5-Gp}-  vpNdv} 


© 

© 


© 


(5.39) 


The  various  contributions  have  again  been  labelled;  a  few  blank  spaces  are  purpose¬ 
ly  indicated  for  comparison  with  the  one-dimensional  results.  Apart  from  the  in¬ 
clusion  of  residual  burning  in  the  chamber,  Eq.  5.39  is  the  same  as  a  result  ob¬ 
tained  in  Ref.  3. 


5  In  this  term  only,  7  •  u  has  been  replaced  by  w _/p,  by  use  of  the  continuity 
in  Eq.  4.1  written  for  the  mean  flow.  This  substitution  aids  comparison  with  the 
one-dimensional  result. 
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These  results  for  the  complex  wave  numbers  will  be  discussed  further  in  later 
sections.  First ,  however,  they  will  be  used  to  motivate  a  fourth  formulation  of 
the  analysis  of  instabilities. 


NWC  1*P  5349 


6.  THE  MODIFIED  THREE-DIMENSIONAL  PROBLEM 


At  the  end  of  Section  4,  it  was  noted  that  the  classical  three-dimensional 
formulation  does  not.  reduce  to  the  one-dimensional  results  covered  in  Sections  2 
and  3.  Moreover,  it  is  not  possible,  within  the  strict  three-dimensional  analysis, 
to  handle  coupling  at  the  boundary  when  the  main  wave  motions  are  parallel  to  the 
boundary-  Loth  of  these  difficulties  are  removed  by  patching  the  three-dimensional 
analysis. 

The  essential  point  is  that  the  one-dimensional  formulation6  implicitly 
accounts  for  viscous  processes  occurring  in  the  flow  adjacent  to  the  lateral  boun¬ 
daries.  If  those  terms  correctly  representing  boundary  layer  effects  can  be  iden¬ 
tified  in  the  one-dimensional  ayalysis,  then  they  may  be  incorporated  in  certain 
results  of  the  three-dimensional  analysis. 

Since,  the  differential  equations  for  the  one-dimensional  problem  contain  these 
boundary  terms  (e.g.,  those  terms  involving  nq,  and  in  Eq.  2. 1-2. 4)  while  the 

partial  differential  equations  for  three-dimensional  problems  do  not,  it  is  neces¬ 
sary  to  establish  a  connection  elsewhere).  In  the  present  work,  it  is  convenient 
to  compare  the  formulas  for  the  complex  wave  numbers  in  Eq.  5.30  and  5.39. 

The  problem  comes  down  to  determining  what  terms  of  Eq.  5.30  represent  genuine 
boundary  effects  and  should  therefore  be  included  as  surface  terms  in  the  three- 
dimensional  result  in  Eq.  5.39.  Resolution  rests  on  correct  interpretation  cf  the 
various  integrals.  Comparison  shows  that  the  one-dimensional  result  in  Eq.  5.30 
contains  six  terras  more  than  the  three-dimensional  result  in  Eq.  5.39.  These  are 
numbered  in  Eq.  5.30  and  will  be  considered  in  turn.  It  must  be  emphasized  that 
the  right-hand  side  of  Eq.  5.30  represents  the  various  contributions  to  rates  of 
energy  addition  to  acoustical  motions  entirely  parallel  to  the  lateral  boundary, 
while  the  mass  addition  (mt,,  £5,  etc.)  is  occurring  mainly  normal  to  the  surface. 
Deviations  at  the  surface  from  motion  normal  to  the  surface  are  represented  by  m 
and  u.  It  should  also  be  noted  that  in  all  terms,  one  factor  p^,  or  in  some  cases 
dp^/dz,  is  present  because  of  the  spatial  structure  of  the  acoustic  mode.  Thus, 
for  example,  a  process  which  couples  to  the  pressure  fluctuation  is  most  effective 
where  p^  is  large.  The  numbered  terms  in  Eq.  5.30  have  the  foil-wing  interpretations. 


e 

The  difference  between  the  formulations  of  Sections  1  and  2  lies  solely  in 
the  way  in  which  variations  of  area  are  handled.  This  is  of  no  consequence  for  the 
present  discussion,  so  the  two  analyses  will  be  referred  to  together  as  the  one¬ 
dimensional  formulation.  It  should,  however,  be  emphasized  that  the  modified  one¬ 
dimensional  formulation  produces  more  accurate  results  when  the  cross-sectional 
area  varies  along  the  axis  of  the  chamber.  The  results  of  the  two  formulations  of 
course  agree  for  the  case  of  constant  area. 
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Q  and@  .  This  represents  the  direct  influence  of  fluctuations  of  the  mass 
flux  from  the  lateral  boundary.  It  arises,  of  course,  due  to  the  coupling  between 
the  acoustic  field  and  the ^combustion  processes  at  the  boundary;  it  should  be  com¬ 
pared  with  the  term  ipak^[up^Sjje]|j  which  represents  the  corresponding  interaction 
at  the  end  surfaces  where  there  is  no  component  of  acoustic  velocity  parallel  to 
the  surface.  In  the  latter  case,  the  driving  of  the  waves  is  proportional  to  the 
velocity  fluctuation. 

The  difference  between  the  two  situations  may  be  interpreted  as  follows. 

When  the  surface  is  oriented  normal  to  the  wave  motion,  motion  of  either  the  sur¬ 
face  or  of  the  gases  leaving  the  surface  produces  p-v  work:  it  is  the  rate  of  p-v 
work,  proportional  to  the  velocity,  which  matters.  But  suppose  that  the  surface 
is  oriented  parallel  to  the  wave  motion.  The  influence  of  fluctuations  of  the 
flow  leaving  the  surface  is  not  now  as  an  agent  doing  work  on  the  waves,  but  as  a 
mechanism  for  locally  increasing  the  pressure  of  the  acoustics  wave.  The  origin 
of  this  term  is  a2  /  mjjdq  on  the  right-hand  side  of  Eq.  3.3.  It  may  be  recalled 
that  Eq.  3.3  is  the  result  of  adding  R/Cv  times  the  energy  equation  (Eq.  3.2)  to 
RT  times  the  continuity  equation  (Eq.  2.1);  that  step  produces  the  source  term 

t^(hos-eoHRT:li'mb<J4  =  CyM»RAT0)J  nyiq  . 

There  are,  in  fact,  two  perturbation  terms  arising  here,  for 

(a2Jmbdq)'  =  a2Jmjjdq+(aZ)'j*inbdq  =  a2  Jm^dq+yRT*  Jir^dq 

=  a2Jmj3dq+a2(^-)  SL  JS^dq  .  (6.1) 

P 

If  only  the  first  term  is  taken  into  account  on  the  right-hand  side  of  Eq.  3.5, 
one  has 

P*SC)+\P-^(U'SC)  =  a2jmjdq  .  (6.2) 


Thus,  the  term  in  question  appears  as  a  source  of  pressure  changes. 


As  an  aside,  it  is  instructive  to  examine  an  approximate  means  of  determin¬ 
ing  the  contribution  from  a  from  Eq.  6.2;  differentiate  Eq.  6.2  with  respect  to 
time  and  replace  3u’/3t  by  the  approximation 


9u* 

at 


JL  9£l 

—  dz 
P 


» 


so 


at 


(p'Sc)-a‘ 


9< Vs  J  . 


8z 


dm! 

=  a2 J  -£jj—  dq  «  ia3k,  J  m^dq 


(6.3) 


Now  suppose  Sc  is  constant  and  note  that  for  standing  waves. 
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£ 


-k5>' 


—2,2  , 

If-  y\* 


O  P‘  — Cy  *•  , 

— C-  «  -a  k  p' 
9t 


so  Eq.  6.3  gives 


(lc2-k^)p 


Now  multiply  by  p'  and  integrate  over  the  chamber  to  find 

^  SJJ  V*-  • 

c  0 

which  is  exactly  that  exhibited  in  Eq.  5.30  for  Sc  constant. 

The  point  is  that  this  term  involving  m^  is  a  source  term  associated  with 
turning  of  the  flow  at  the  boundary;  it  should  therefore  also  appear  in  the  three- 
dimensional  case.  However,  it  should  be  present  only  when  there  is  a  component 
of  the  acoustic  velocity  parallel  to  the  surface.  To  see  what  should  be  used  in 
the  general  three-dimensional  problem,  it  is  helpful  to  contrast  the  two  cases 
which  have  arisen  in  the  one-dimensional  problem.  The  contributions  to  (k2-k2)E|, 
per  unit  area,  for  the  direct  coupling  between  the  acoustics  field  and  combustion 


no  acoustic  velocity 
parallel  to  the 
surface 


.  A  A  .  /(,  — 2 

•ip  akjPjU-i  — -  upj 

a 


(6.  5a)7 


acoustic  velocity 
only  parallel  to 
the  surface 


-iaV A-1  -  RhA  K  *  (6* 5b) 


It  should  be  recalled  that  the  acoustic  velocity  referred  to  here  is  that  of 
the  classical  unperturbed  acoustic  modes.  At  the  surface  itself,  the  normal  com¬ 
ponent  is  always  zero  for  those  motions;  the  u  in  Eq.  6.5a  and  nq>  in  Eq.  6.5b  are 
the  small  fluctuations  normal  to  the  surface  due  to  coupling  between  the  acoustic 
field  and  combustion.  Term  ©  has  also  been  included  in  Eq.  6.5b  since  it  is 
also,  obviously,  associated  with  coupling  at  the  boundary;  it  too  comes  from  the 
source  term  (a2  +  RAT0)  /  m^dq. 

7  The  minus  sign  arises  when  the  first  term  in  Eq.  5.30  is  written  out;  e.g., 
[up^Sije]|}  =  (SbePjju)^  -  (SjjePu)0-  At  z  =  0,  u  is  positive  for  flow  outward  from 
the  surface,  which  is  the  convention  followed  In  Eq.  6.5a  and  6.5b. 
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The  question  which  must  be  faced  here  concerns  the  correct  form  of  the  coupl¬ 
ing  terms  when  the  unperturbed  acoustic  velocity  in  the  vicinity  of  the  surface  is 
neither  purely  normal  to  the  surface  [case  (Eq.  6.5a)]  nor  purely  parallel  to  the 
surface  [case  (Eq.  6.5b)].  Several  different  arguments  have  been  examined,  but  it 
appears  that  the  desired  result  can  best  be  obtained  in  the  following  way.  By 
definition,  the  fluctuation  of  mass  flux,  m  =  pu,  is 


~  +  £- 
Uu  Pv 


“b 


(6.6) 


so  the  velocity  fluctuation  is 

m 


£  =  — =  ~  (IVnV?r)  * 

“b  p  p  p 


(6.7) 


By  use  of  the  perfect  gas  law,  the  density  fluctuation — which  here  is  at  the  edge 
of  the  combustion  zone — is 

o  _£  T 

~  _  -r-  •  (6.8) 

P  P  T 

The  temperature  fluctuation  can  be  broken,  as  earlier,  into  two  parts:  that  asso¬ 
ciated  with  the  wave  motion  outside  the  combustion  zone — the  isentropic  fluctua¬ 
tion — and  the  nonisentropic  part: 


T  T  is. 

Thus,  the  density  fluctuation  is 


T_  _  /  T  \  +  AT  _  y-1  p  +  AT 

rVi  m  ' *  'm  "V  r  '•Ti 


(6.9) 


and  Eq.  6.7  becomes 


A  _  A  .  _ 

_£_  =  J_  _  AT. 

P  V  ?  T 


1  -  -  1  P  ,  *%. 

u  =  ~  mb“Ub7i‘+=? 

p  P  P T 


AT  . 


(6.10) 


(6.11) 


The  last  result  is,  of  course,  always  true,  irrespective  of  the  character  of  the 
acoustic  mode.  Now  substitute  this  into  Eq.  6.5a  to  find 

. - %  -  Ub*,mbA*>1 

-ipak.p. - p .  + - AT  1  -  i  —  up.  . 

^  P  VP  PT  "  a  4 
The  convection  terms  cancel,  and  the  result  is 
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no  acoustic  velocity 
parallel  to  the 
surface 


l 


■iak^«.mb-^~  RptATmb 


(6. 12) 


which  should  be  compared  with  Eq.  6.5b.  Thus,  the  only  difference  between  the  two 
cases  is  the  factor  y  in  the  nonisentropic  term  in  Eq.  6.12. 


Now  for  the  general  case,  the  acoustic  velocity  of  the  unperturbed  acoustic 
field  near  the  surface  is  partly  parallel  and  partly  normal.  The  problem  of 
representing  the  surface  response  in  this  case  is  discussed  in  Appendix  A,  where 
it  is  shown  that  Eq.  6.5b  and  6.12  should  be  weighted  respectively  by  factors 
g||,  6X,  such  that  6||  +  6x  =  1;  in  the  limit  of  no  parallel  velocity,  6||  =  0, 

6x  -  1»  while  if  the  acoustic  velocity  is  entirely  parallel,  6||  =  1,  6x  =  0. 

These  weighting  factors  are  related  to  the  mode  structure  as  shown  in  Appendix  A. 
The  discussion  there  breaks  essentially  into  two  pieces:  a  detailed  examination  of 
a  general  acoustic  mode  in  the  vicinity  of  a  rigid  surface;  and  a  "proof"  for  the 
representation  of  the  coupling  in  the  general  case  (Eq.  A-15). 

The  conclusion  of  the  above  argument  is  that  the  missing  terms  and  Q 
in  Eq.  5.39  are  accounted  for  and  combine  with  the  surface  combustion  terms  al¬ 
ready  present,  in  the  surface  integral 

Rm, 

_iakN  §PN[«ib  +  ~-z —  Af(6||+y8x)]dS  .  (o.l3) 

a 

The  remaining  boundary  terms  in  Eq.  5.30  will  be  incorporated  in  the  three- 
dimensional  problem  in  a  similar  way. 


2  .  This  term  represents  an  energy  exchange  associated  with  the  acquisition, 
by  the  average  inflow,  of  the  local  acoustical  energy.  Again  it  appears  as  a 
source  term  in  the  one-dimensional  problem  and  is  present  only  if  there  is  a  com¬ 
ponent  of  acoustic  velocity  parallel  to  the  surface.  Hence,  the  corresponding 
term  for  the  three-dimensional  problem  is 

E$n  +“T  <vPn)2^  6>i“bdS  *  (6*14) 

pa  kjj 
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Q)-  This  term  represents  an  energy  loss  associated .with  turning  of  the 
average  flow  of  particles  into  the  chamber  from  the  boundary.  The  corresponding 
term  in  three  dimensions  is  evidently 


5  UP  dz 


J  mjp^dqdz 


'  ^X^I'dS  . 


(6.15) 


0-  This  term  involves  fluctuations  of  the  inflow  or  outflow  itself,  paral¬ 
lel  to  the  boundary.  It  can  therefore  be  taken  over  directly  to  the  three- 
dimensional  problem.  At  the  burning  surface,  both  p  and  p  may  reasonably  be 
taken  to  be  zero — i.e.,  the  flow  is  assumed  to  depart  the  surface  in  the  normal 
direction. 


All  of  the  preceding  arguments  lead  to 

<k2-kN,EN  *  {-iIkN$pN£’V'Sr  (6|i+y6x)  ws) 

a 

+  {i^f  Cpn+Jz  (Ppi/lSniVs-i^  JpN2”“dv} 

pa  kN  ^  a 

+  {§up  .  Vp^H^dS-  JF  •'C'p^v} 

+  {-i  —  -3-  J  PNGe_w +e _w  +(Q+Q  )]dV 
L  a  CV  N  P  P  P  P  P 


J?w  dV+  J 
N  p  J 


w 


A  A  v  —  A 

(S-S 

>  P 


N 


dV 


} 


-  §[^111^  +  up,|mb(p)]-  7pNdS  .  (6.16) 

This  is  a  new  result;  in  later  sections  it  will  be  applied  to  some  specific 
problems.  It  should  be  noted  first  that  because  of  the  way  in  which  Eq.  6.16  has 
been  constructed,  it  reduces  exactly  to  the  ■modified  one-dimensional  result  in 
Eq.  5.30  Also,  it  should  be  recalled  that  some  of  the  boundary  terms  have  been 
combined  and  written  in  a  different  way.  The  modified  one-dimensional  result 
corresponding  exactly  to  Eq.  6.16— i.e.,  if  Eq.  6.2  is  used  in  place  of  Eq.  6.5  is. 


(k2.kt2)E^  = 


* 

AT 


P<Xe- 


L 

0 
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"iak^  0  ^1‘‘fI?5bdq+:^>fA^mbdq]dS} 
+ v  =■  J*  frt +tt  Car )  ]  J™bdqdz 


k.  L  o  , 

,P  1  H 

_i~  j  p*  th  (^sc)scdz 

a  0  &c  az  c  c 


+  {“i  —  J*P.  [e  w  +e  w  +(Q+Q  )]S  dz+i J p?w  S  dz 

L  T  Cvn  '  P  P  P  P  P  c  ::vt  pc 


L.dP, 


a  VO 


o  *  p 


L  dp ,  •.  L  dp . 

+  J  (u -u  }w  —r-^  S  dzl  -  J  £  S  dz  . 
g  ppdz  cJg^dz  c 


(6.17) 


In  summary,  it  is  perhaps  helpful  to  emphasize  what  has  been  accomplished. 
What  is  eventually  required  is  the  formula  for  the  wave  number  of  three- 
dimensional  oscillations.  But  the  purely  three-dimensional  formulation,  for 
inviscid  flow,  fails  to  contain  certain  contributions  arising  essentially  from 
the  flow  in  a  boundary  layer.  Strictly,  those  contributions  can  be  found  only  by 
actually  solving  the  boundary  layer  problems.  However,  the  one-dimensional  form¬ 
ulation  affords  a  means  of  determining  approximately  the  desired  results.  The-"----* 
accuracy  of  the  results  cannot  be  established  using  the  procedure  followed  here. 

The  idea  is  usefully  illustrated  by  a  more  familiar  example,  the  attenuation 
of  waves  due  to  viscous  forces  in  the  acoustic  boundary  layer  on  the  wall  of  a 
tube.  This  is  usually  found  by  computing  the  energy  dissipation  in  the  (linear¬ 
ized)  acoustic  boundary  layer.  For  a  standing  wave  in  a  circular  tube  of  diameter 
D,  the  contribution  to  a  is 

a  =  -  ^  Vv~i  ( 1  -f  Irl—  f  (6.18) 

D  7pF 

where  v  is  the  kinematic  viscosity  and  Pr  is  the  Prandtl  number.  This  result  is 
based  on  the  assumption  that  the  boundary  layer  is  locally  one-dimensional,  with 
the  velocity  at  the  outer  edge  equal  to  the  total  acoustic  velocity. 

On  the  other  hand,  there  is  a  displacement  effect  of  the  boundary  layer, 
which  can  be  taken  as  the  basis  of  a  different  calculation  which  has  been  worked 
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out  by  Chester.  Because  the  fluctuation  of  velocity  varies  along  the  axis  of  the  ' 
tube,  3u/3z  ^  0,  preservation  of  continuity  in  the  boundary  layer  requires  that 
the  velocity  v,  normal  to  the  axis,  at  the  edge  of  the  boundary  layer,  be  non-zero. 
If  the  acoustic  equations  are  written  in  one-dimensional  form  for  the  cross-section 
of  the  tube  bounded  by  the  edge  of  the  boundary  layer,  then  one  has  Eq.  2. 1-2. 4 
except  that  there  are  now  no  particles,  no  combustion  of  any  sort,  and  the  cross- 
sectional  area  is  fixed.  The  only  difference  from  the  usual  one-dimensional  equa¬ 
tions  is  a  source  term  in  the  continuity  equation. 


|^  +  ~  (pu)  =  sLJpvbdq  =  ^-  pv 


(6.19) 


There  is,  of  course,  no  average  flow  inward,  so  vb  has  only  a  fluctuating  part. 
Chester  has  computed  vb  by  analyzing  the  boundary  layer  flow,  and  the  linearized 
form  of  Eq.  6.19  is 


■ar  +  p Is-  =  P  =  w  (6‘20) 


*/Pt  0 


vs 


while  the  momentum  and  energy  equations  are  simply 

—  3u‘  3p|  _ 

p-sr+it  =  0 


(6.21) 


-  3T’  p  3u'  „ 

P-5T+C  -5T  =  ° 

V 

It  follows  that  the  wave  equation  for  the  pressure  fluctuation  is 

,2  .  . 


-El  _  JL  _  1 

£  -2  Tr  ~  "v 


1  9w 


9z**  a°  91“  ^ 


(6.22) 


(6.23) 


For  harmonic  motions  having  complex  wave  number  k,  the  procedure  used  in 
Section  5  then  produces  the  formula  for  a  standing  wave  in  a  tube. 


2  2  im  0 

k  =k/„“-r 


L 

J  wp  dz 


,2  .  2®  ^  . 
=  k2 ' 1  yTj  J  w?< 


dz 


(6.24) 


The  imaginary  part  gives  the  formula  for  o, 

a  =  vL  l  w  P>iz 


(6.25) 
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Now  for  a  standing  wave. 


ik.  iu).t 

- -  e  sin(k.z)  . 

pa 


9n 

0z 


iU)^t 

e  cos(k  z) 

<v 


and 


id) 


w  =  we 


' =  ^  C1 + ^)(=?)cos(kii,e 


CO 


r  e-iU)$ 
o  yf 


of  which  the  real  part  of  w  is 


z)  . 


Since  p^  =  cos(k£z),  Eq.  6.18  gives 

a  =  -  -$■  a/vuTT  (l  +  }  .  (6.26) 

D  1  V  VP?  ' 

This  differs  from  Eq.  6.18  only  in  the  replacement  of  by  2/y.  For  air 
(y  =  1.4)  the  error  is  very  small  indeed — much  less  than  can  be  measured,  in 
fact. 


Although  a  boundary  layer  analysis  was  used  to  evaluate  vj,  in  Eq.  6.19> 
this  example  does  show  how  a  one- dimensional  formulation  with  sources  at  the 
lateral  boundary  may  reproduce  a  result  found  wholly  by  a  boundary  layer  analysis, 
i.e.,  that  give  Eq.  6,18. 
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7.  A  FEW  REMARKS  ON  INTERPRETATION 


In  Section  5,  Eq.  5.6,  the  growth  constant  a  was  related  to  the  rate  of 
change  of  the  time-averaged  total  energy  in  the  chamber.  It  is  helpful  to  inter¬ 
pret  the  formulas  in  Section  6  in  that  way.  The  real  and  imaginary  parts  of  Eq. 
6.16  are 


-kjX  =  {ak  §  PNCA«>  Af  W(8|+y«xn«} 

a  a 

+{§  “pr’-',%6«SbP,dS-I ^<r)'  7PNdV} 

+{rr  #•  f  Ph^VVp111*®  (i,+Q^)]dv 


_  ^  *  (r)  -> 

+Jw  (u-u  )  *  vPNdv} 


-  fS  c  ^ii(r^b + ^r,^p)  3  *  ^PN^S 


(7.1) 


-2a  ^  4  =  -{^N  &'  PNCAbr>  +  Af  (r)(Siiky6J.)]ds} 

a  a 


+  {~Z  §  tPN+  TT  (^PN^lSilir^dS  -  ~  J  P^7*  udv} 

pa  kN  a 

+  {  #'u(i).V^6n^dS  -  J*  F(i)*  7^dv} 
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+  ^-JPN*'pdV  +  J'"p("'“o,(  •VPNdV} 


The  fact  that  ujj  is  pure  imaginary  (see  Eq.  5.36)  has  been  used  in  Eq. 
7.2. 


(7.2) 

7.1  and 


It  is  mainly  Eq.  7.2  which  is  of  practical  interest,  and  especially  the 
first  five  terms  representing  the  coupling  between  wave  motions  and  combustion, 
and  the  mean  flow/ acoustics  interactions.  In  order  to  clarify  the  meaning  of 
these  terms,  it  is  helpful  to  express  o  as  the  ratio  of  the  rate  of  change  of 
acoustic  energy  to  the  total  acoustic  energy  in  the  chamber,  as  in  Eq.  5.6. 

First  note  that  E^  is  proportional  to  the  total  time-averaged  acoustic  energy  in 
the  chamber.  Thew instantaneous  time-averaged  energy  density,  with  the  time 
factor  dropped,  is: 


e,T  = 


'N 


*  2 
1  ,  PN 
4 

pa 


2 

+  P^n  )  =  p 


N 


+  K. 


N 


(7.3) 


where  PN  and  %  are  the  time-averaged  densities  of  potential  and  kinetic  energy. 
Hence,  the  total  time-averaged  energy  in  the  chamber  is 


dV 


(7.4) 


Now  by  use  of  the  wave  equation  for  pjj, 

J>PN)2dV  =  JP-(PNPPNMV-J'PN72%dV 


=  SdS  +  kN2  JPNdV 


The  first  integral  vanishes  according  to  the  boundary  condition  on  pN,  so  that  the 
normalization  integral  is  related  to  the  time-averaged  total  energy  in  the 
chamber  according  to 


2pl2e. 


N 


(7.5) 
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TKe  formulas  for  the  one-dimensional  formulation  are  entirely  analogous,  giving, 
for  example. 


5jtScdz  -  EU  -  2Pa  ^ 


(7.6) 


with  the  time-averaged  energy  density 


-  2 

,Po  _  .  2. 

‘imk  (w** 


(7.7) 


By  taking  the  imaginary  part  of  Eq.  6.17,  one  finds  for  the  growth  constant 
in  the  one-dimensional  problem 


afc<iiA”br,+''  A,?(r)pJ  -|X[J“br,d<l+:§  J'AT(r,mbdq]d2i} 

2r  LI  3 


a  p  0 


I.1-  2  1  d  No  . 

’f  P*  s“  dz  (uSc)Scdz 
0  c 


+ K  t ^upr>VV'pr)+<Q<r>+Qpr,)]scd2 
+ v  {  n"vscAz+ I’-pC5'1  )-iprl)  ~a£  SM 

■  -J  Wi^Jmhdq+SpW  J^'dq]  -|§  dz 


(7.8) 
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To  see  what  the  various  pieces  mean,  it  is  best  to  replace  E§  by  use  of  Eq.  7.4, 
Ele  ky  Eq.  7.6,  and  use  the  fact  that  -2a  is  the  fractional  rate  of  change  of  time- 
averaged  energy  in  the  chamber  (cf.,  Eq.  5.6): 


2a  =  111  •  (7*9) 

For  -2a§N  and  -2aC£,  one  finds  the  sum  of  terms  shown  in  the  following  table. 

TABLE  1. 

-2a£.  I  -2aG T 
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TABLE  1.  (Cont'd.) 


©  M 


P.w  S  dz 
l  p  c 


-= -  fwD(S(r,-S<r>  &  S  dz 

2p  ak,  0  P  P  dz  c 


© 

& 


2-y  f  P,Tw  dV 
1  J  N  p 


r~  J ■w„(“<r,-uir,  )-7pn«v 


2PakN 

-1 


2p^kN 


r~*  (i)—  -*  (i)—  (p)-,  * 

Fuh  p>vPNds 


A  factor  ^  appears  explicitly  in  terms  (2)  (4)  (5)  (8)  and  (ll)  because  of 

the  time  averaging.  For  example,  consider  thg  first  term  in  ©■  If  Eq.  6.12  is 
used  to  eliminate  mgr'  in  favor  of  and  AT(r),  a  term 

i<r>  Sbe!C 

a  a  (r) 

arises.  The  product  p^u  is  the  rate  at  which  work  is  done  per  unit  area  of 
surface,  and  the  h  is  the  average  of  cos  ut  cos  (ait  +  <j>)  where  <f>  is  the  phase 
difference  between  p^  and  u'  at  the  surface.  The  remaining  pieces  of  the  groups 
(D  and (2)  are  also  rates  of  doing  work,  associated  essentially  with  fluctuations 
of  density.  Similarly,  all  other  terms  containing  the  \  are  various  kinds  of 
time-averaged  rates  of  doing  work. 

The  term  ©  shows  that  there  is  a  time-averaged  energy  loss  at  the  surface, 
equal  (per  unit  area)  to 

-  2  -  2 

y  (PN+KN^b6|i  "  2  (~3-  +  IP  ~f")ab6l! 

The  factor  6 II  ,  as  argued  above,  assures  that  this  term  arises  only  if  there  are 
acoustical  motions  parallel  to  the  surface.  As  emphasized  already,  the  process 
represented  by  this  term  is  the  inelastic  acceleration  of  the  gases  leaving  the 
surface  normally,  so  that  they  acquire  the  acoustical  motions  parallel  to  the 
surface.  The  time-averaged  energy  possessed  by  the  gas  is  ultimately  (Pjj  +  Kjj). 
The  factor  2  arises  here  because  the  amount  of  work  done  in  such  a  process  is 
twice  the  mechanical  energy  finally  gained.  Half  of  the  work  done  appears  as  the 
mechanical  energy,  and  half  is  dissipated  as  heat.  The  total  work  is  of  course 
done  by  the  waves  already  in  the  chamber  and  hence  is  obviously  an  energy  loss. 

Similarly,  the  term  (l^)  represents  an  inelastic  process;  the  gases  evolved 
in  the  combustion  of  particulate  matter,  at  the  rate  Wp,  must  acquire  the  local 
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potential  energy,  Pjj.  Once  again  the  work  done  is  twice  the  energy  finally  pos¬ 
sessed  by  the  gas. 

Finally,  term  ©  should  be  interpreted  as  a  time-averaged  work  per  unit 
volume  and  time*  For  the  acoustic  field,  pN/p  =  (pN/p),  so  p^/pi2  =  pN/p. 

Hence,  the  factor  in  parenthesis  is  just  Pjj(^*u)/p  and  represents  a  rate  of  change 
of  volume  due  to  the  density  fluctuation  and  dilatation  by  the  mean  flow.  Thus, 
the  whole  term  does  indeed  represent  a  rate  of  work  by  the  pressure  fluctuation. 

In  summary,  it  is  possible  to  provide  a  physical  interpretation  for  all  the 
terms  in  the  growth  constant,  although  their  appearance  is  not  obvious  in  all 
cases. 
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8.  AN  ALTERNATE  METHOD  OF  COMPUTING 
THE  GROWTH  RATE 


In  Ref.  2,  Jantrell  and  Hart  have  proposed  a  means  of  computing  the  growth 
rate  for  acoustic  waves  in  a  rocket  chamber.  Part  of  their  results  are  identical 
with  part  of  the  results  found  here,  but  their  formula  for  the  growth  constant  a 
is  incomplete.  Moreover,  if  their  procedure  is  followed  exactly,  it  is  not  possi¬ 
ble  to  determine,  among  others,  the  influence  of  residual  combustion.  Their  tech¬ 
nique  is  therefore  not  applicable  to  liquid  rockets,  for  example,  whereas  the 
present  analysis  is  valid  for  any  combustion  chamber. 

It  is  the  purpose  of  this  section  to  summarize  the  calculation  proposed  in 
Ref.  2;  to  identify  the  respects  in  which  it  differs  from  the  present  work;  and 
finally  to  show  how  the  general  idea  proposed  there  can,  when  correctly  developed; 
produce  results  in  accord  with  those  found  here.  However,  it  should  be  emphasized 
that  since  the  analysis  of  Ref.  3  is  for  an  inviscid,  three-dimensional  flow  field, 
it  does  not  and  cannot  contain  the  boundary  processes,  discussed  in  Section  6, 
which  have  been  incorporated  in  the  modified  three-dimensional  problem. 


8.1.  SUMMARY  OF  THE  ANALYSIS  OF  REFER.  CE  2 


The  main  idea  proposed  by  Cantrell  and  Hart  is  that  the  stability  of  waves  in 
a  chamber  can  be  studied  by  examining  the  time  rate  of  change  of  energy  in  the 
chamber.  This  is  computed  by  integrating  the  energy  equation  over  the  volume  of 
the  chamber,  and  making  use  also  of  the  equations  of  conservation  of  mass  and 
momentum.  Residual  combustion  and  particulate  matter  were  not  accounted  for  in 
Ref.  2,  but  they  will  be  included  here.  Then  the  inviscid  equations  of  motion 
are  Eq.  4.1,  4.4,  and  4.5.  In  order  to  clarify  the  correspondence  between  this 
work  and  Ref.  2,  Eq.  4.4  and  4.5  are  slightly  modified  to  introduce  the  enthalpy 
h  =  /  CpdT: 


||  +  V  •  (m)  =  wp  (8.1) 

ff  +  V(h  +  =  £  (F  -  a)  (8.2) 

2  2 

h  p(cvT  +  r>  -  -v-iS(h  +  r»  +  vHo"p  +  (9+V  (8-3) 
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where  m  =  pu  is  the  vector  mass  flux  of  gases.  Equations  8. 1-8. 3  are  equations 
(2)- (4)  of  Ref.  2  except  for  the  additional  terms  due  to  residual  combustion  and 
the  presence  of  particles. 

Now  the  left-hand  side  of  Eq.  8.3  is  the  time  rate  of  change  of  the  total 
energy  density  for  the  flow  field:  it  includes  not  only  the  energy  of  acoustic 
waves  but  also  the  energy  associated  with  the  mean  flow.  As  noted  below,  this  is 
a  crucial  observation.  Thus,  integration  of  Eq.  8.3  over  the  mean  volume  of  the 
chamber  gives  the  time  rate  of  change  of  total  energy  withi..  the  chamber: 


2  2 

/  p(CvT  +  ~~)dV  =  ~/V‘[m(h  +  j-)]dV  +  /  [up*F  +  epQWp  +  (Q+Qp)dV. 


(8.4) 


It  is  necessary  to  extract  from  this  equation  an  equation  for  the  balance  of 
acoustic  energy  only. 

To  do  this,  the  variables  are  as  usual  expressed  as  sums  of  mean  values  and 
fluctuations.  Cantrell  and  Hart  show  that  the  left-hand  side  and  the  first  vol¬ 
ume  integral  on  the  right-hand  side  can  be  rewritten  and  partially  combined  to 
yield  their  equation  (16).  Let  <  >  denote  time  averages,  and  one  has  from  Eq. 
8.4  with  their  result  written  using  the  notation  defined  here: 


2p  a  a 


-  .2 


(  $§  dSn*  [ u'p'  +  — --g-  +  p(u*  u')u'  +  (u*  u' )u  ]  ) 

pa  a 


+  <  J  Cup-  F+  epowpo  +  (Q+QpndV  > 


(8.5) 


There  are  now  two  important  points  to  be  made:  (1)  the  last  term  on  the  left- 
hand  side  of  Eq.  8.5  and  the  last  two  terms  in  the  surface  integral  on  the  right- 
hand  side  should  be  dropped  under  the  circumstances  for  which  the  analysis  is 
useful;  (2)  the  last  set  of  terms  in  Eq.  8.5  will  not  produce  the  influences  of 
residual  combustion  and  particulate  matter  found  in  the  present  work,  for  a 
reason  given  below. 

As  Cantrell  and  Hart  have  argued,  the  left-hand  side  might  be  interpreted 
as  the  expression  for  the  time-averaged  rate  of  change  of  total  acoustic  energy 
in  the  chamber.  Hence,  without  the  differentiation,  it  is  the  time-averaged 
acoustic  energy,  and  the  right-hand  side  represents  the  time-averaged  rate  of 
change  of  energy  due  to  the  mean  flow/acoustics  interactions,  residual  combustion, 
and  particulate  matter.  Therefore,  by  the  definition  in  Eq.  7.9  of  the  growth 
constant  a. 
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-2a  = 


—  ,  2  ,  -* 

{  £j>  dSn*  Cp'u'  +^2"  +  p^u*  tf  )u'  +  — (u»  u*)u  ]  ) 

_ e_s _ z± _ 

2p  a  a 

(JdV[Sp-Ftepowp-MQ4-Qpn> 

<Jdv[ffS'2  +  ^+-S#)P']> 

2p  a  a 


(8.6) 


Unlike  the  calculations  given  in  the  preceding  sections  of  this  report,  Eq. 
8.6  is  formally  valid  to  any  order  of  the  mean  flow  Mach  number.  Consequently, 
p*  and  u'  represent  acoustic  quantities  which  also  may  differ  from  the  classical 
unperturbed  values  by  terms  of  first  and  higher  order  in  the  Mach  number.  How¬ 
ever,  for  practical  purposes,  it  is  certainly  a  convenience,  quite  likely  ade¬ 
quate,  and  virtually  a  necessity  if  relatively  simple  results  are  obtained,  to 
consider  perturbations  only  to  first  order  in  the  Mach  number. 

At  the  burning  surface,  the  first  approximation  to  u'  (which  is  of  course 
zero  there  for  the  unperturbed  field)  is  of  the  order  of  the  mean  flow  Mach 
number,  as  for  example  shown  in  the  definitions  in  Eq.  9.1  and  9.2  for  pressure 
coupling;  |A^|  and  | R^, J  are  around  1  -  7  or  so.  If  the  amplitude  of  th|  oscilla¬ 
tion  is  taken  t<j  be  of  order^unity^then  p'u'  is  of  order_TL,  as  is  up’  /p32  in 
Eq.  8.6.  But  (u'u^u1  and  (u  •  u’)u  are  no  greater  than  M^.  Hence,  if  those 
terms  are  retained,  higher  order  terms  in  p’  and  u'  must  also  be  retained.  In 
order  to  find  those,  one  must  return  to  the  inhomogeneous  differential  equations 
and  solve  them,  a  step  which  has  never  been  taken  in  practice,  and  obviously  must 
involve  a  great  deal  of  effort.  (Formulas  for  the  first  correction  to  the  un¬ 
perturbed  values  have  been  given  in  Ref.  3.)  For  the  same  reason,  the  term 
(u  •  u^)p’/ai  should  be  dropped  from  the  denominator  of  Eq.  9.6,  since  it  is  of 
order  Mj,  and  therefore  adds  a  correction  term  to  a  of  order  Mb2.  Thus,  the  first 
part  of  Eq.  9.6  is  correctly 


The  denominator  is  precisely  Efj2/2pa2  and  this  result  is  identical  to  that  found 
from  the  real  part  of  the  first  term  in  Eq.  5.39. 

It  is  obviously  desirable  to  avoid  computing  the  acoustic  field  to  higher 
orders;  numerical  procedures  are  available  for  computing  the  classical  field  in 
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any  geometry,  but  not  with  perturbations  due  to  the  mean  flow.  This  is  the  main 
reason  why  regions  of  high  Mach  number,  such  as  exhaust  nozzles  and  the  vent  of  a 
T-burner,  should  be  treated  separately;  they  are  then  introduced  in  the  acoustics 
problem  as  boundary  conditions,  represented  by  admittance  functions. 

Now  consider  the  second  part  of  Eq.  8.6.  It  is  necessary  to  expand  the  quan¬ 
tities  to  second  order,  in  whatever  small  parameters  are  present;  for  example. 


=  F  +  + 


The  first  ordj^r  tejpns  are  associated  with  the  unperturbed  acoustic  field.  Consid¬ 
er  the  term  Up  •  F 


(u  +u 
P  pl 


.  )•  (f+f1,+f2,+. 


) 


=  u  *F+(u  *Fi  +u'  ,*  F )+(u 


1 


Pl 


F'+u'  ->• 

1  p2 


F+u 


The  first  term  is  associated  with  the  average  flow  and  should  be  dropped.  The 
time  average  of  the  first  order  terms  is  zero.  In  this  way,  the  denominator  of 
the  ratio  gives 


<  I  dVC<V  FI,  J+(ep0iWpj+ep02Wp+Sp0Wp2^(Q2+Qp2 )]  >  <8'  8> 


Not  only  are  the  contributions  from  particulate  matter  in  disagreement  with  those 
found  here,  and  used  elsewhere,  but  the  influence  of  heat  release  appears  only  in 
second  order.  The  second  result  is  clearly  wrong;  for  example,  such  terms  must 
arise  to  first  order  if  the  behavior  of  liquid  and  gas  rockets  is  to  be  described 
correctly.  The  first  order  effect  does  appear  in  the  result  in  Eq.  6.16  found 
here. 


The  question,  then,  is:  why  does  the  analysis  of  Ref.  2  produce  the  correct 
form  for  the  mean  flow/acoustics  interactions  (apart  from  the  boundary  terms) 
but  fail  to  represent  correctly  the  influence  of  particulate  matter  and  residual 
combustion?  The  answer  follows  from  the  fact  that  the  analysis  is  based  on  Eq. 
8.4  for  the  time  rate-of-change  of  the  total  energy  in  the  chamber,  rather  thin 
on  the  energy  balance  for  the  acoustic  field  alone,  although  the  latter  is 
apparently  obtained  eventually. 

For  by  using  the  total  energy,  one  cannot  in  any  obvious  way  show  exchanges 
of  energy  occurring  totally  within  the  system.  Thus,  for  example,  heat  lost 
from  the  mean  flow  by  combustion  and  added  to  the  acoustic  waves  appears  twice; 
once  with  a  +  and  once  with  a  -,  thus  cancelling.  The  debit  and  credit  accounts 
are  not  shown  in  detail,  but  only  their  net  result.  A  net  result  does  appear  at 
the  boundary,  as  shown  by  the  first  term  of  Eq.  9.6. 
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8.2.  THE  EQUATION  OF  BALANCE  FOR  THE  TOTAL 
ENERGY  OF  THE  ACOUSTIC  FIELD 


The  idea  of  Ref.  2,  that  the  growth  constant  can  be  found  from  considerations 
of  the  total  acoustic  energy,  can  be  developed  to  produce  results  in  exact  agree¬ 
ment  with  those  found  here.  The  starting  point  is  the  linearized  equations  in 
Eq.  4 . lO^and  4.11  for  the  three-dimensional  acoustic  field.  Take  the  scalar  pro¬ 
duct  of  u'  with  Eq.  4.10  and  add  the  result  to  p'/pl2  times  Eq.  4.11  to  find 


J^fp  " =T"  *  2  )  +  (U'P'^  =  -P^'*  ( a#  Vu1 +u'*  Vu )  +  -  %  g  u- Vp'+u'*  (F'-a'J+p'P* 

'  2p  a  pa 

P  (8.9) 


Now  integrate  over  the  mean  volume  of  the  chamber  and  convert  the  integr  '1  of 
V*(u'p')  to  a  surface  integral: 


d£ 

dt 


+  dSn*  (u'p')  =  -f  dVrp  u1*  (u*  Vu'tu*-  Vu)  + — ^  u»  Vp1  ]  +JdVu'-  (F'-a1 ) 


P  a 


JdVp'P' 


(8.  10) 


This,  of  course,  is  valid  only  to  first  order  in  the  average  Mach  number,  and  one 
can  use  unperturbed  values  for  u'  and  p'  on  the  right-hand  side. 

The  terms  in  the  first  integral  on  the  right-hand  side  can  be  combined  in 
essentially  the  same  way  as  those  following  Eq.  5.38  were  treated.  By  use  of 
some  vector  identities. 


u' •  (u*  Vu'+u'- Vu)  =  u' •  [V(u- u' )-u'XVXu' ]  =  u'.V(u-u') 

=  V*  L(u»  u' )u'  1  -  (u-  ul )V-  u1  . 

Thus,  the  volume  integral  is 

-4  — *  ^  -4 

JdV[  pV*  {{uvu1  )u' }  —  p  (u*  u')V*  u'+ ***  VP*  1 

P  a 

-4  «-4  j2  “♦ 

=  dSn*  u'p  u*  u'+JdV[-p  u»  (u'V-  u)+ u*  Vp*  1 

P  a 

~  O 

The  integrand  of  the  surface  integral  is  of  order  and  hence  must  be  dropped. 
For  harmonic  motions,  the  first  term  in  the  second  integral  can  be  written  by 
use  of  the  acoustic  equations 

“♦  1  -4 

u'  =  — _ ■  Vp'  ,  V*  u 

ip  ak 


iak 

YP 


P' 
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so 

-7u*  (u'V-u')  =  u-  Vp')  =  ^2  u*  7p' 

\p  P a 

which  equals  the  second  term  in  the  second  integral,  as  demonstrated  in  the  argu¬ 
ment  following  Eq.  5.38.  Thus,  the  volume  integral  reduces  to 


— *  — *  ^  1 2 

JdV[  p  u'»  ( u«  Vu'-Hi1-  Vu)-f  u»  Vp1  ]  =  .fdVu.  V 

pa  pa 


pa  pa 


«2  r 


=  0ds5-u(^)-Jdv^v.u 

pa  pa 


which  is  the  result  found  just  preceding  Eq.  5.39. 

Now  to  obtain  this  result  in  general,  it  is  necessary  to  show  that 

|  — * 

-J* dVfTu*  (u* V.  u' )  =  JdV  u  *  Vp* 

pa 

or 

J  dV u*  [  Vp*  +  u’V.  u'  ]  =  0  . 
pa 


Again  by  use  of  the  classical  acoustics  equations,  this  can  be  written 

JdVu-  [P'-~£  -*’-^3  =  o 

and  by  differentiating  with  respect  to  time,8 


at 


-  a  V 
a? 


8  Because  of  this  step,  the  end  result  in  Eq.  8.11  has  been  shown  to  be  true 
only  to  within  a  constant  which  can  clearly  be  taken  to  be  zero  with  no  loss  of 
generality. 
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But  both  u*  and  p'  satisfy  the  unperturbed  wave  equation,  since  terms  of  order 
must  be  neglected;  hence 


3^u’  3^p'  _  | — 2-*  — 2  .  n 

p'  — y  -  u'- — t-  =  -p  a  u'  +u'a  p'  =  0  . 

at  at 

Therefore,  the  desired  combination  is  obtained,  i.e., 

-*  — *  f  -*  ”•  ,2 

JdVt-'pu*  (u'V.  u'  VP'3  =  JdV  Z —Z  u*  Vpl  =  JdVu*  ——t  ) 


P  a 


pa 


-»  ,2  ,2  -» 

dSn*  u(  P ~2  )  -  JdV  u  (&• 

pa  pa 


and  Eq.  8.10  can  be  put  in  the  form 

2  2  -» 

^  dSn-  [u'p' +u-£-^-  3+JdVu'- (F’-a'J-JdV;^- V*  u+JdVp'P' . 

pa  pa 

(8.  12) 

The  corresponding  value  of  a  is  given  by  the  formula  obtained  after  taking 
time  averages  of  Eq.  8.12 

-  ,2  ,2 

-2cr£  =  <  <$)  dSn  -  [u'p'+  i^  ])  -  (  JdV  -E~  V-  u  >  +  (JdVu1’  (F'-a')> 

pa  pa 


+  (  JdVp'P'  > 


(8.  13) 


It  is  a  simple  matter  to  verify  that  for  harmonic  motions  Eq.  8.13  is  identical 
with  the  real  part  of  Eq.  5.39,  which  of  course  was  obtained  by  an  entirely 
different  computation. 

Hence,  it  has  been  shown  that  by  considering  the  balance  of  acoustic  energy 
one  can  obtain  the  same  result  for  the  growth  constant  as  that  found  by  an  approx¬ 
imate  solution  to  the  linearized  equations  of  motion.  However,  it  is  important 
to  note  that  both  analyses  cannot  provide  the  surface  terms,  essentially  associa¬ 
ted  with  viscous  boundary  layer  effects,  which  have  been  incorporated  in  Section 
6  in  the  modified  three-dimensional  analysis. 
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9.  A  SIMPLE  APPLICATION  TO  T-3URNERS 


It  is  an  important  result  of  the  analysis  developed  here  that  not  only  is 
the  response  of  a  burning  surface  characterized  by  two  functions,  but  both  neces¬ 
sarily  appear  in  a  stability  analysis.  Consider  the  case  of  pure  pressure  coupl¬ 
ing;  then  the  admittance  (A^)  ard  response  (R^)  functions  are  defined  by  the 
formulas 


a  VP 


A 


m  VP 


(9.1) 


(9.2) 


These  relationships  apply,  of  course,  at  the  edge  of  the  combustion  zone.  The 
signs  are  set  so  that  if  A},  and  Rj,  are  real  and  positive,  u  and  m^  are  real  and 
positive  outward  from  the  surface.  After  substitution  of  these  definitions  into 
Eq.  6.11,  one  finds  an  expression  for  AT  in  terms  of  Aj,  and  R^: 

AT  r^H)  ..  TJ  1 JL 

___  =  r— -  +  l-RbJ-^  .  (9.3) 

T  Mb  VP 


The  combination  in  brackets  might  also  be  called  the  entropy  response  func¬ 
tion  for  pressure  coupling,  since  At/T  =  As/C„  where  As  is  the  amplitude  of  en¬ 
tropy  fluctuations  at  the  edge  of  the  combustion  zone.  To  see  this,  begin  with 
the  combined  expressions  of  the  First  and  Second  Laws  of  Thermodynamics 

Tds  =  de  +  pdv  =  de  -  — 

P  P 

3ut  for  a  perfect  gas,  de  =  CvdT  and  dp/p  =  dp/p  -  dT/T,  so  one  has 

Tds  =  C  dT-^  =  C  dT-RT^  . 

P  P  P  P 

Now  apply  this  to  small  fluctuations  of  an  element  of  gas,  and  split  dT  into  a 
sum  of  isentropic  and  nonisentropic  parts: 

dT  _  +  AT 

T  VP  T 
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Substitution  in  the  previous  equation  leads  to 


Tds  =  CpAT  +  Cp-R]|  dp  =  CAT 


(9.4) 


Since  (y-l)/yCp  =  R,  the  second  term  vanishes,  leaving  the  result  quoted. 


The  implications  of  these  results  may  be  most  easily  seen  by  examining 
the  simple  case  of  a  T-burner  having  samples  of  propellant  at  the  ends,  normal  to 
the  axis,  and  short  grains  al$o  at  the  ends,  but  extending  along  the  lateral 
walls.  It  will  be  assumed  that  the  latter  are  approximately  flush  with  the 
chamber  wall  so  that  the  cross-sectional  area  of  the  chamber,  Sc,  is  constant. 
Moreover,  resi’ual  combustion  and  the  presence  of  particulate  matter  will  be 
ignored  so  that  Eq.  7.8  simplifies  to 

.2a^i  =  ^  Af'r,pt)Sbef 

a  a  0 


-q  J ^br )+  ^  AT(r  ,mb]dz} 

0  a 

k-t5  p  r.2  1  ,dPt  fn-  ktsc  Kzdn, 

—  J  [p<.+7^  (ir  >  ]mbdz  -  —  L  pt  dSdz 


(9.5) 


a  0 


The  perimeter  of  the  chamber  is  q,  so  that  if  the  length  of  the  lateral  grains  is 
Lb,  the  area  of  the  burning  surface  on  the  sidewalls  is  S]jS  =  qL},. 

With  the  definitions  in  Eq.  8.1  and  8.2,  the  combustion  terms  in  Eq-  8.5  are: 

Rni  jj 

a-ZktSbe(A^,+SV  <9-6) 

a  u 

and 

-plk^qj  . 

For  short  grains,  p^  -  1  under  the  integral,  and  with  the  assumption  that  the  re¬ 
sponse  functions  are  constant  over  all  the  burning  surface,  the  second  contribu¬ 
tion  is  simply 


(9.7) 


-2ktsbSrR^r,+  7  <^^)+Mb-MbRl^r, «  • 
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The  factor  2  arises  in  both  Eq.  9.5  and  9.6  because  grains  are  placed  at  both  ends 
of  the  chamber;  the  burning  area  of  one  lateral  grain  is  S^s  =  qLfc,  where  is 
the  length  of  the  grain. 

Since  dp^/dz  =  0  near  the  ends,  the  surface  interaction  terms,  evaluated  for 
the  lateral  burning  surfaces  only,  are 


pa 


Wbs 

"pa 


2k, S  M 
l  bs  b 


(9.8) 


This  acoustics/mean  flow  interaction  term  also  contributes  at  the  vent.  For  the 
fundamental  (and  all  odd  modes),  -  0  there,  while  dpj/dz  =  -k£  sin(k£L/2)  = 
-k£.  Hence,  at  the  vent,  x 


k.  L 


p  a  0  k^ 


A 

v 


where  \  is  the  flow  area  of  the  vent,  and  Sy  is  the  magnitude  of  the  mass  flux 
through  the  entrance  of  the  vent.  The  sign  change  occurs  because  the  flow  is 
outward.  Now  by  continuity. 


mvAv  "  2mb(Sbs+Sbe’ 
so 

“-2WSbe+Sbs'  •  «•» 

pa  0  k^  vent 

This  holds  for  all  modes  in  the  simple  uniform  burner,  since  for  the  even  modes, 
p^  -  1  at  the  center  of  the  burner,  but  dp^/dz  ~  0,  and  one  finds  again  Eq.  9.9. 

Note  that  there  is  a  gain  of  energy  associated  with  the  flow  out  the  vent. 

It  arises  because  as  the  flow  exhausts,  there  is  an  inelastic  process  at  the 
boundary  in  which  the  fluid  originally  participating  in  the  acoustical  motions 
flows  out  the  vent  with  no  acoustic  energy.  The  energy  it  loses  is  given  up  to 
the  field  in  the  chamber,  thereby  constituting  a  gain.  This  process  is  distinct 
from  the  radiation  of  acoustic  energy,  a  loss  which  is  associated  with  pressure 
fluctuations  at  the  vent,  and  is  therefore  present  only  for  the  odd  modes.  It 
should  be  emphasized  that  this  contribution  from  the  vent  is  a  formal  result  of 
the  one-dimensional  analysis  and  in  no  way  constitutes  a  detailed  analysis  of 
the  vent.  Indeed,  one  should  anticipate  that  the  one-dimensional  approximation 
is  much  less  accurate  in  the  vicinity  of  the  vent,  so  that  the  results  are  sus¬ 
pect.  Moreover,  while  it  is  relatively  easy  to  understand  that  there  should  be 
an  energy  loss  associated  with  the  influx  of  mass  at  the  surface,  it  is  difficult 
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to  understand  how  the  inverse  process  actually  takes  place  at  the  vent.  This  ques¬ 
tion  is  an  important  one  for  the  interpretation  of  data  taken  from  T-burner  fir¬ 
ings,  but  it  will  not  be  examined  here  in  detail.  A  careful  consideration  of  the 
vent  should  probably  produce  a  response  or  admittance  function  which  appears  as  an 
additional  contribution  to  the  energy  balance  and  hence  the  growth  constant. 

The  last  term  in  Eq.  9.5  can  be  rewritten  by  using  the  continuity  equation 
for  the  mean  flow, 

psc  pS5 


Over  the  burning  surfaces  one  finds 


k  S 
-f  c 

a 


A. 
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aZ  du 

pl  Hz 


dz 


?  Vbs”^ 

pa 


2k  ,S,  M 
l  bs 


b 


(9.10) 


This  term  also  provides  a  contribution  at  the  vent,  but  for  the  even  modes  only. 
Since  p^,2  =  1  there. 


?+_  u 

1  du  s»2k.S  —  . 

n  l  C  — 

0-  a 


Here  u  is  the  axial  velocity  in  the  chamber,  and  again  by  continuity, 

Sc"  =  <Sbs  +  Sbe^b 


so  that 


r  k  S  L  o  j —  1 

a  0 


dz  1  «  2k^(Sbe  +  Sbs)Mb 


(9.11) 


for  even  modes  only. 

Finally,  for  this  special  case, 

2 


L  o 

U)  it  p  *  2„  .  JLj  ,  _ 

=2EU  =  -=  I  h  ScdZ  “  ~2  Vc  • 

a  a  u 


(9.12) 


Substitution  of  Eq.  9.6-9.12  into  Eq.  9.5,  and  some  rearrangement,  gives 
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%■  --  {2(A^r,+Mb,T£  +2tMbRbr,+7l |H 


K<¥^>-K'' 


S,  +S  -i  CL-L 

Hb,-^,]  }.-A- 

c  even  a 

modes 
only 


(9.13) 


The  term  a<jL/a  represents  other  possible  sources  of  attenuation,  such  as  wall  fric¬ 
tion.  The  first  set  of  brackets  contains  all  contributions  from  the  coupling  of 
pressure  fluctuations  with  combustion;  the  second  set  contains  the  mean  flow/ 
acoustics  interactions  at  the  burning  surface  and  at  the  vent. 

An  interesting  special  case  is  that  for  odd  modes  (really  only  the  fundamen¬ 
tal  is  of  practical  interest);  if  the  coupling  should  happen  to  be  isentropic,  then 
MbRb  =  Ab+Hb  and  Eq.  9.13  becomes 


3  L  /_v  S,  +S,  a,l 

—  =  2(Aj|r)+2Mb)(-^|-  --)  --1- 
a  °c  a 


(9.14) 


The  2  multiplying  Mb  is  half  due  to  the  end  discs  and  half  due  to  the  vent.  This 
result  has  motivated  testing  with  the  variable  area  T-burner.  Measurements  of  the 
growth  constant,  a,  are  made  for  different  area  ratios,  (Sbe+Sbs)/Sc.  Then  if  a 
is  plotted  versus  the  area  ratio,  the  slope  is  2(ab^+2Mb),  and  the  intercept  is 
ad*  This  procedure  is  based  on  the  crucial  assumption  that  the  attenuation  aj  is 
independent  of  area  ratio,  a  circumstance  which  is  likely  not  to  be  true. 

More  importantly,  however,  there  is  no  justification  for  assuming  that  the 
coupling  is  isentropic.  But.  if  that  assumption  fails,  Eq.  9.13  contains  essen¬ 
tially  three  unknowns:  Ab^r',  Rb^ri  and  otjj.  Again  for  odd  modes,  one  has 


2(A^rW 


b'Ly  S, 


+S,  a  .L 

+21^  (-bg-—  )  -  -4 


(9.15) 


In  principle,  a  sequence  of  three  tests  (at  the  same  frequency,  say)  in  which 
Sbe  an<*  ^bs  are  var*e(*  will  provide  three  values  of  a,  and  hence  Eq.  9.15  leads 
to  three  equations  in  the  three  unknowns  A^uO+Mb,  Rb^r\  and  a<j.  It  remains  to 
be  seen  whether  such  a  strategy  may  be  useful. 

Real  T-burners  in  which  the  grains  are  not  flush  with  the  lateral  walls  re¬ 
quire  more  elaborate  analysis  outside  the  scope  of  this  work. 
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10.  APPLICATION  OF  THE  MODIFIED  THREE-DIMENSIONAL 
RESULT  TO  MODES  IN  A  CYLINDRICAL  CHAMBER  AND 
THE  DATA  OF  BROWNLEE  AND  MARBLE 


The  most  extensive  study  of  a  stability  boundary  in  solid  propellant  rockets 
is  that  due  to  Brownlee  and  Marble  (Ref.  5).  In  Ref.  6  and  later  in  Ref.  3,  the 
classical  linearized  stability  analysis  was  applied.  Although  attenuation  by 
particles  was  suggested  as  the  dominant  damping  mechanism  in  Ref.  6,  viscous 
stresses  at  the  head  end  were  supposed  to  be  more  important  in  the  discussion  of 
Ref.  3.  However,  not  until  much  later  was  the  admittance  function  of  the  propel¬ 
lant  measured  in  a  T-bumer  (Ref.  7). 

Midway  in  the  experiments  program  reported  in  Ref.  7,  it  appeared  that  the 
analysis  of  Ref.  3  had  been  verified  (Ref.  3).  Subsequently,  this  turned  out  not 
to  be  the  case.  The  propellant  used  in  the  experiments  of  Ref.  5  contains  no 
metal,  so  it  is  unlikely  that  attenuation  by  particles  in  the  product  gases  is 
important.  Indeed,  motion  pictures  show  that  the  products  are  in  fact  very  clean. 
There  appear  now  to  be  at  least  two  possible  reasons  why  the  analysis  and  measure¬ 
ments  do  not  agree — the  discrepancy  is  quite  large  (Ref.  7): 

(1)  The  T-burner  measurements  of  Ref.  7  provide  only  values  of  the 
admittance  function; 

(2)  The  various  boundary  terms  introduced  in  the  modified  three- 
dimensional  analysis  have  until  now  not  been  taken  into 
account. 

Unfortunately,  there  is  still  no  data  for  the  response  function.  Hence,  the  ques¬ 
tion  of  the  nonisentropic  behavior  associated  with  the  difference  between  the  re¬ 
sponse  and  admittance  function  cannot  be  explored.  The  purpose  here  is  simply  to 
examine  the  importance  of  the  new  terms  introduced  in  the  modified  three-dimen¬ 
sional  analysis. 

All  influences  of  condensed  particles  and  residual  combustion  will  be  ig¬ 
nored,  but  the  possibility  of  nonisentropic  processes  at  the  boundary  will  be  in¬ 
cluded.  Only  the  stability  of  waves  will  be  treated  here,  so  the  problem  is 
simply  to  evaluate  the  right-hand  side  of  Eq.  7.2: 
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2  *=T  En  =  {a  $  Af  ^r^(6|ffv6x)]ds] 

^  ct 

-{=:  $  tp^+-z  (^leiiSvds} 
pa  kN 

+  iT'  S“b“»<i)- '?PNds - 2  ^  en  •  uo.i) 

N  a 

As  earlier,  all  other  processes  of  attenuation  have  been  represented  by  a  decay 
constant,  a<j. 

In  the  results  reported  by  Brownlee  and  Marble,  a  particular  mode  was  ob¬ 
served,  but  there  is  no  added  work  here  in  computing  the  growth  constant  for  all 
modes  in  a  cylindrical  chamber  and  later  specializing  the  results  as  required. 

The  classical  acoustic  modes  are  given  by  the  formula: 

pN  =•  cos(kjz)  cos(m4>)  Jm(Kmnr)  (10.2) 


^-=H 

P  akjq  P  a*cN 


UN  ‘  ~  rN 


dJ  A 

r— —  cos(k,z)cos(m4>)-4>“~  J  cos(k  z)sin(m<j>) 

QX  «v  X  IH  *v 


-z  k^  Jmsin(k^z  )cos 


(10.3) 


2 

+  x 
mn 


(10.4) 


In  order  to  satisfy  the  boundary  condition  on  pjq,  xmn  are  the  roots  of 

d  J  ,  , 

(x  r)  =  0  (r  =  r  ) 

dr  mn  c 


where  rc  is  the  radius  of  the  chamber.  As  usual,  m  =  0,  1,  2,  ...  and  k£  =  £v/L 
with  £  =  0,  1,  2  ...  The  mean  velocity  field  for  incompressible  potential  flow 
is  _ 

u  =  u,[  2  —  z  -  r  1  •  (10.5) 

"  C 
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It  is  simply  a  matter  of  arithmetic  to  compute  V*um  and  hence 

*  tv* 

7 |j  .  Uj^.  ,  6j_ ,  and  6j|.  For  the  cylindrical  chamber. 


«  XNT 

V±’  ^  3r 

*  f  -  1  9UN<I>  ,  aUNs 

"‘UN  r  3r  3  z 


(10.6) 


(10.7) 


In  these  computations,  some  use  is  made  of  the  equation  satisfied  by  the  Bessel 
function, 

d2J  ,  dj  ?  2 

- y—  + - j —  +  (k  - )J  =  0  .  r io.  8) 

,2  r  dr  mn  2  m. 
dr  r 

At  the  surface  (r  =  rc)  one  then  finds 


vr  uT 


*  2  2  2 

-  -  — z  -Wc  )JmC0S<V)c0s(m'l', 


(10.9) 


-  2  2  2 

VU-  u.  =  -  (m  - k^  r  )Jmcos(kjZ)cos(m<J>) 

<-Vc  '  ‘ 

and  therefore,  by  substitution  into  Eq.  6.19  and  6.20, 

,2  2  2, 

m  -k  r 

c  „  1  mn  c  1  I 

bx  ‘  l  "j  z  2 — r. - r~T~rr  J 

1x11  -  Wc  l  +  'm  -Vc  I 


(10.10) 


(10.11) 


2.22 
m  -k,  r 

•t  c 


,  2  2  2,'  ,  2,22 

m  +  m  -k.  r 

'  nan  c  1  '  t  c 


1 

2T  J 


(10.12) 


These  formulas  show,  as  one  would  anticipate  from  the  symmetry  of  the  situation, 
that  the  weighting  factors  5fl  and  Sx  are  independent  of  position  on  the  surface. 
Note  that  only  for  purely  radial  modes  (m  =  =  0)  does  8||  vanish. 


With  these  preliminary  calculations,  it  is  not  difficult  to  evaluate  the 
contributions  to  a. 
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10.1.  COUPLING  AT  THE  BURNING  SURFACE. 


It  will  be  assumed  that  only  pressure  coupling  is  important.  With  AT  given 
by  Eq.  9.3,  the  <:  upling  terms  are 

-  (r  j 

^[Rbr  ^6ll+6^>+7  (~-  +1-Rbr>)]  JL{Hmnrc)*[COs2(in*)d<t>|Cos2(k^z)dz 


YP 


■  MbSbs[Rir)<5»«^7  (^- 


(r) 


(10.13) 


where  S^  =  2-nr^L  is  the  area  of  the  burning  surface,  and 


W  =  ?  <1+W1+6*o>  • 


(10.14) 


The  Kronecker  deltas,  6mo,  6^0,  are  unity  if  m  =  0,  Z  =  0,  respectively  and  zero 
otherwise. 

10.2.  ACOUSTICS/MEAN  FLOW  INTERACTIONS  AT  THE  BURNING  SURFACE 

This  term  is  non-zero  only  along  the  burning  surface.  Although  6 1|  is  non-zero 
at  the  entrance  plane  of  the  nozzle,  and  of  course  there  is  a  mean  flow  there, 
there  is  no  contribution.  The  reason  is  that  this  term,  by  construction,  is 
present  only  if  the  acoustic  field  is  zero^  outside  the  boundary  considered;  for 
only  then  must  there  be  a  transition  through  a  boundary  layer.  The  energy  loss  is 


—  JJ  CPn+~T^Pn*  ^5||mbdS  ” 

P a  kjj 

,  r  mZ+kirZ  i 

=  Mbsbs5" V Vcl1  + — zrr^JAmi  .  (10-15) 

rc  *N 


10.3.  INFLUENCE  OF  THE  NOZZLE 

The  first  Surface  integral  of  Eq.  10.1  also  applies  at  the  entrance  to  the 
nozzle.  Relatively  little  is  known  about  the  nozzle  losses,  but  formally  at 
least,  admittance  and  response  functions  can  be  introduced  in  the  same  way  as 

^  Or,  more  generally,  discontinuous.  But  non-zero  values  just  outside  the 
boundary  different  from  those  just  inside  will  not  occur  in  practice. 
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for  the  burning  surface.  The  values  of  6j.  and  6||are  of  course  different;  for  ex¬ 
ample,  Vj_.  is  now  3ujjz/3z.  One  then  finds 


l-  l  „2  i  ,,  2  2 

2m  -x  r  +k  r 
mn  c  1  -t  c 


(10.16) 


2  2  2 1 
2m  -  k  r  I 
1  mn  c  1 


mn  c 


2  2  2.2 
2m  -k  r  +k,  r 
1  mn  c  1  l  c 


(10.17) 


The  signs  on  the  mass  fluxes  are  changed,  ffij,  -*■  -me  and  mj,  -*■  me  because  the  flow 
is  outward;  then  the  response  function  for  the  nozzle  is,  for  example. 


_£  =  R  JB.  . 
—  n  — 

a*  yp 


(10.18) 


The  contribution  of  the  nozzle  is  then 


b  0 

It  is  likely  that  the  admittance  function  (in  particular  the  real  part)  is  very 
small  unless  the  mode  has  an  axial  component  of  velocity.  If,  in  addition,  the 
coupling  is  isentropic,  Eq.  10.19  reduces  simply  to 


<1+6mo>J  „r)rdr 


mo  g  m  mn 


(10.20) 


which  represents  a  loss  of  wave  energy  due  to  convection  by  the  mean  flow  exhaust¬ 
ing  from  the  chamber. 

The  next  to  last  term  in  Eq.  10.1,  depending  on  u  contributes  nothing. 

At  the  burning  surface  it  is  zero  because  the  flow  leaving  the  surface  may  be 
assumed  to  have  initially  no  motion  parallel  to  the  surface.  At  the  nozzle  en¬ 
trance,  there  is  no  contribution  because  the  parallel  motions  (i.e.,  normal  to 
the  axis)  are  continuous  through  the  entrance  plane. 


Finally,  the  values  of  Ejj  are  given  by 

3*  ^ 

E*  =  2ttLA  .  J*  J2  (x  r)rdr  =  S  r 
N  m-t  ^  m  mn  °s  c  c 


-^--Ar2 

2  2/  n 
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Bringing  all  of  the  above  results  together,  one  has  for  a: 

E2  r  a(*) 
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(10.22) 


(other  losses) 

Consider  the  mode  reportedly  observed  by  Brownlee  and  Marble,  namely  that 
for  m  ®  1  and  Z  =  n  =  0.  For  the  nozzle,  6n|j  =  0  and  6^  =  1 »  while  along  the 
burning  surface  both  6||  and  6x  are  non-zero.  The  Mach  number  at  the  nozzle 
entrance  is  found  from  Eq.  10.5  to  be  He  =  2 (L/rc)b^ .  Eventually,  Eq.  10.22 
can  be  reduced  to 

*olrc 


(x2r2+l)  A(r) 

-6,|  -|4—  .  (_S_  +  i)(. 


1+8 


to 


)- 


a,r 
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(10.23) 
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On  the  stability  boundary,  a  =  0,  and  Eq.  10.23  can  be  rewritten  as 

(H2,r2+1K 

c  ol  c  '  \ 

-6 1|  5  2 - I  • 


(10.24) 


‘■Vi1;-*) 


For  isentropic  flow,  Rb  =  Ab  /Mb  +  1  ,  and  if  both  the  nozzle  admittance 
function  and  the  loss  associated  with  mass  addition  at  the  burning  surface  are 
neglected,  the  Eq.  10.24  reduces  to  Eq.  24  of  Ref.  6;  for  Jt  =  0, 


(r) 


a.r 
d  c 


Ml  x  ,r  aM, 

ol  c 


1 

~z — Z 
*01*0 


(10.25) 


Now  \>lrc  =  »  so  the  second  term  on  the  right-hand  side  of  Eq.  10.25 

is  -.295.  As  argued  above,  and  in  Ref.  6,  losses  within  the  volume  are  apparent¬ 
ly  negligible;  but  that  leaves  only  viscous  damping  on  the  head  end  as  the  source 
of  loss.  If  one  assumes  simple  laminar  shear  losses,  it  turns  out  that  the 
first  term  in  Eq.  10.25  is  about  one-tenth  that  of  the  second.  Hence,  specifica¬ 
tion  of  the  stability  boundary  comes  down  to 


=  -  .295  ,  (10.26) 

M. 

D 

which  represents  a  simple  balance  between  acoustics/combustion  coupling  and 

a c oust ic 53 /mean  flow  interactions.  Measured  values  of  Ab  ' /Mb  ,  taken  from 

T-burner  data  by  using  Eq.  8.13  for  end  discs  only  (S^e  =  Sc)  are  shown  in  Fig. 

10.1  taken  from  Ref.  7,  for  p  =  300  psia. 

Note  that  Ab^r^  <  0  for  frequencies  greater  than  approximately  2000  Hz,  which 
means  that  the  combustion  processes  are  absorbing  energy.  The  -.295  is  the  net 
driving  effect  by  the  mean  flow,  shown  as  the  dashed  line  In  Fig.  10.1. 

By  taking  data  for  several  mean  pressures,  a  plot  of  the  intersections  of 
Ab^/Mb  with  -.295  gives,  according  to  Eq.  10.26,  the  stability  boundary  shown 
in  Fig.  10.2.  There  is  clearly  a  very  large  difference  between  this  result  and 
the  observations  of  Brownlee  and  Marble. 

Now  suppose  that,  while  An  ^  is  still  zero,  the  coupling  is  not  isentropic, 
and  the  loss  due  to  turning  of  the  flow  is  taken  into  account.  Define  AR^  as 
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ARb  = 


Then  Eq.  10.24  leads  to 


(10.27) 


^-  +  (^i)6,|ARb 


k  ,r 
ol  c 


Holrc 


(10.28) 


If  ARb  £  0,  this  may  provide  a  different  result  for  the.  "predicted”  sta  ility 
boundary.  To  the  present  time,  measurements  of  Rb  (i.e.,  T-burner  tests  with 
grains  on  the  lateral  boundary)  have  not  been  taken,  so  it  is  not  possible  to 
determine  whether  or  not  the  explanation  proposed  here  is  adequate  to  account 
for  the  poor  agreement  shown  in  Fig.  10.2.  Indeed,  another  very  strong  possibil¬ 
ity  is  that  the  response  of  the  propellant  to  the  azimuthal  acoustic  velocity 
parallel  to  the  burning  surface — i.e.,  velocity  coupling — may  be  significant. 
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APPENDIX  A 


The  problem  discussed  here  sesms  to  be  one  which  has  not  previously  been 
treated.  It  is  concerned  with  the  question  of  representing  in  general  the  coupling 
between  a  surface  and  an  ur-'teady  flow,  accounting  for  the  fact  that  the  represen¬ 
tation  depends  on  the  flow-tield  itself.  From  the  analysis  of  the  stability  prob¬ 
lems,  one  finds  that  the  interaction  between  an  acoustic  field  and  a  burning  sur¬ 
face  depends  cn  the  component  of  velocity  parallel  to  the  surface.  When  Eq.  6.5b 
and  6.12  are  inserted  in  the  expressions  for  the  complex  wavenumber,  one  finds 


2  ,  2.  2  . 
{k  =  < 


-iakt 

a 


ail  parallel  (A- la) 

acoustic  velocity 


no  parallel 
acoustic  velocity 


(A- lb) 


-iak^  J]  frn^  +y  ^  n^AT  3  p^dS 
a 

when  only  these  coupling  terms  are  considered. 

To  motivate  the  argument,  take  the  imaginary  part  of  Eq.  A-l  to  find  for  the 
growth  constant 


El  P 

«"> 

(A-2a) 

^  §  C-3  I^r)+ VltAT (r)3p^dS 

(A-2b) 

V.  P. 

The  integrands  have  the  units  of  velocity;  write  them  respectively  as 
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a  (r) 

U|| 

+  ^RAT  {r) 

P 

(A- 3a) 

a  (r) 
uj. 

+  u^yRAT  (r) 

P 

(A- 3b) 

2  -2 

As  shown  in  Sectlor.  7,  Eg  =  2pa  £■£ 
energy  in  the  chamber;  hence,  Eq.  A 

where  is  the  time-averaged  total 

-2a, b  can  be  written 

acoustic 

7E"  JT  p^ui/r)dS 

(A-4a) 

2a  =  ^ 

Trt  XT  p^x(r  >dS 

(A-4b) 

Obviously,  then,  the  t  .me-averaged  rate  of  work  done  on  the  waves  by  the  inter¬ 
actions  at  an  element  of  surface  dS  is 


(A-5a) 


(A- 5b) 


It  should  be  emphasized  that  for  each  case,  the  fluctuations  u||  ,  ux  are 

normal  to  the  surface;  the  subscripts  denote  only  that  feature,  of  the  unper¬ 
turbed  acoustic  field,  used  to  distinguish  the  two  limiting  cases. 

That  these  different  results  have  been  obtained  for  the  coupling  in  the  two 
limiting  cases  implies  that  the  coupling  for  the  general  case  somehow  depends  on 
the  relative  orientatj-m  of  the  surface  element  and  the  local  acoustic  field,  for 
this  determines  what  Lhe  parallel  velocity  will  be  for  a  given  pressure  amplitude. 
To  develop  this  idea,  observe  first  that  however  complicated  the  mode  structure 
may  be,  the  field  in  the  vicinity  of  a  surface  may  be  approximated  locally  by  a 
superposition  of  incident  and  reflected  waves.  Let  9  be  the  angle  between  the 
normal  to  the  surface  and  the  directions  of  those  waves.  It  is,  of  course,  the 
unperturbed  acoustic  field  near  a  rigid  surface  which  is  of  concern  here.  Sche¬ 
matically,  the  limiting  cases  and  an  arbitrary  acoustic  mode  may  be  represented 
as  shown  in  Fig.  A-l.  For  all  parallel  velocity,  9  =  rr/ 2  and  for  no  parallel 
velocity,  €  =  0.  Note  that  for  the  unperturbed  field,  the  normal  component 
u  •  n  of  the  velocity  is  zero. 


(rate  of  work)  = 


1  £  „  (r )  iq 

2 


1  *  ,*  (t)jQ 

2  P,ul  aS 
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ALL  PARALLEL 

ACOUSTIC 

VELOCITY 


NO  PARALLEL 

ACOUSTIC 

VELOCITY 


WAVE  TO  RIGHT 


WAVE  TO  LEFT 


INCIDENT  WAVE 


REFLECTED  WAVE 


GENERAL 

ACOUSTIC 

MODE 


WAVE  INCIDENT 
FROM  LEFT 


WAVE  INCIDENT 
FROM  RIGHT 


FIG-  A-l.  Sketch  of  Acoustic  Modes 
in  the  Vicinity  of  a  Rigid  Surface. 


Let  n  and  s  denote  coordinates  normal  and  parallel  to  the  surface  as  shown. 
Then  for  the  wave  incident  from  the  left,  the  incident  and  reflected  waves  may  be 
represented  as 


Pt'  '  =  -pacp’(s  sin  9-n  cos  ft-  at) 

=  -  sin  ftcp’(s  sin  9-n  cos  0  -  at) 


u.  '  =  cos  0cp’(s  sin  9-n  cos  9  -at) 
in 


(A- 6) 


=  -  p  afp' (s  sin  9  +  n  cos  9  -  at) 

=  -  sin  9 cp’ (s  sin  8  +  n  cos  0  -  at) 


=  -  cos  0cp'(s  sin  8  4-  n  cos  0  -  at) 


(A-7) 
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Here,  subscript  i  denotes  incident,  r  denotes  reflected,  and  n,  s  denote  the  norma 
and  tangential  components  of  the  velocity.  There  are  corresponding  formulas  for 
the  wave  incident  from  the  right: 

-p  a  \Jf * (— s  sin  0  -  n  cos  0  -  at) 
sin0  ik*  (-s  sin  0  -  n  cos  0  -  at) 

cos  0  41  (-s  sin  0  -  n  cos  0  -  at)  (A- 8) 


-p  a  ty‘(-s  sin  0  +  n  cos  0  -  at) 
sin  0  i!r'(-s  sin  0  +  n  cos  0  -  at) 

-  cos  0  il'’(-s  sin  0  +  n  cos  0  -  at)  (A- 9) 


The  functions  <P  and  ?  are  velocity  potentials;  the  primes  denote  differentiation 
with  respect  to  the  argument:  e.g.,  cp*  =  dfp(z)/dz.  (See,  for  example.  Ref.  9, 
p.  266.)  Note  that  Eq.  A*-6,  A-7,  A-8  and  A-9  represent  solutions  to  the  unper¬ 
turbed  wave  equation,  satisfying  the  condition  that  the  velocities  normal  to  the 

surface,  u(^)  +  u(£)  and  u(r)  +  u(r)  vanish  at  the  surface  (n  =  0). 
in  rn  in  rn 

Now  expand  u^  =  u^P  +  u^jp  in  Taylor  series  about  the  origin,  along  the 
n  axis,  and  retain  only  the  first  term: 

=  cos0LrP(-ncos0-at)-<p(ncos0-at)3 

s=0 

as  cos0n  [cp(-ncos0-at)-cp(ncos0-at)]n_Q  ~  -2cos  ncp'(-at)  . 


/u(^  =  (u.W+u 

\  n  /  \  in  rn  / 


Similarly,  for  the  wave  incident  from  the  right. 


-2  cos 


2 


On*' (-at)  . 


Thus,  the  sum  of  these  gives 


2Q  1  /^n\  _ 

cos  6  S3  -  \n/n^0  ~  "  2(cp'+^) 


(A-10) 
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Here,  un/n  can  be  interpreted,  since  un  =  0  at  n  =  0,  as  VA*u,  the  divergence  of 
the  velocity  taken  in  the  direction  normal  to  the  surface.  To  see  this,  note 
that  u  =  un  +  us  and  V  operates  only  on  un. 

Similarly,  expand  the  velocity  parallel  to  the  surface,  u|^  =  +  u^g^, 

and  u(r)  to  find 
s 

(us)  "  W  **  -2sin2s(cp'+i|r*)  . 

n=0  ' s=0 

n=0 


In  the  limit  s  *>  0,  [ns(s)-us(0)] /s  can  be  interpreted  as  V||-u,  so 

-  2q  _  vr u 

sm  0  "  ‘  2(qp'+V) 

It  follows  from  Eq.  A-10  and  A-ll  that 


tanB  = 


(A-ll) 


(A- 12) 


where  absolute  values  have  been  taken  to  assure,  as  necessary,  that  only  positive 
numbers  arise.  Then  Eq.  A-12  implies 


*  i 


JT>\\  =  sinQ  =  | 

|V||*u|  J 

_  A  A  J 

|V|i-u|  +  j V^  u  j 

»/5 1  =  cosfl  = 

[-  1 

|V||*u|  +  j  Vi*  u  | 

(A-  i3) 


(A- 14) 


It  may  be  noted  that  Eq.  A-13  and  A-14  are  likely  not  to  be  the  most  useful 
forms  for  practical  geometries,  since  generally  the  pressure  field  is  computed 
numerically.  However,  by  using  the  momentum  equation,  Vp'  =  -p  3u'/3t,  one  can 
write,  for  example. 


Vi*  u 


p(6n)-p(0) 

(bn)2 


V||.  u 


1 

iptu 


p(s+6s)-p(s) 

(8s)2 
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where  6n,  6s  are  the  incremental  distances  measured  in  the  normal  and  tangential 
directions  from  the  point  n  =  0,s.  In  this  way,  the  results  of  a  numerical  calcu¬ 
lation  can  be  used  directly  to  compute  6y  and  6^. 

But  now  the  question  to  answer  is,  "how  should  the  limiting  cases  shown  in 
Eq.  A-l  be  combined  to  give  the  correct  representation  for  the  coupling  when  0  is 
neither  0  nor  tt/2?"  An  appealing  guess  is  that  the  coupling  for  no  parallel  velo¬ 
city  should  be  multiplied  by  cosz0,  and  that  for  all  parallel  velocity  bysinz0. 
That  this  is  correct  may  be  demonstrated  as  follows.  First,  owing  to  the  way  in 
which  the  one-dimensional  analysis  of  Sections  2  and  3  is  formulated,  it  is  clear 
that  it  is  the  component  of  velocity  normal  to  the  velocity  in  the.  chamber  which 
matters.  For  0  #  0,  the  proportion  sin  0  of  the  mats  flux  leaving  the  surface  is 
normal  to  the  acoustic  velocity,  and  the  proportion  cos  0  is  parallel  to  tne  acous¬ 
tic  velocity,  as  shown  in  Fig.  A— 2.  Hence,  a  factor  sin  0  should  multiply  the 
fluctuation  u  (of  which  the  real  part  appears  in  Eq.  A-3a)  or  Eq.  6.5b,  for  tfce 
coupling  due  with  flow  parallel  to  the  surface.  And,  a  factor  cos  0  must  multiply 
U|| ,  or  Eq.  6.12. 


FIG.  A-2.  Projections  of  Velocity  Fluctuations  Associated 
With  Mass  Addition,  and  an  Element  of  Surface  Area. 


But  in  addition  to  weighting  the  velocity  components,  the  surface  areas  must 
be  weighted.  For  the  work  done  due  to  coupling  with  the  parallel  flow,  only  the 
area  dS  sin  0  parallel  to  the  direction  of  the  incident  waves  is  effective,  as 
indicated  in  the  above  sketch.  And  for  the  normal  flow,  the  projected  area  dS 
cos  0  normal  to  the  direction  of  propagation  is  involved.  Consequently,  the  total 
rate  at  which  work  is  done  is 

p(u|x  ^  sin^fi  +  ccs^8)dS 
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It  follows  that  the  two  cases  shown  in  Fig.  A-l  should  be  combined  as 

(jc2_lc2)E  ^=  -iaic  JJ  [m.  (cos20+sin29)+^  m^ATtsin  0+ycos  9)ldS 

X .  Xj  'v  »» 


'l  ** 

=  -ilk.  jJtirih+^-mhAf  (6||+v6x)]dS  (A- 15) 


2  2 

with  6 1.  =  sin  9,  6^  =  cos  6. 


It  must  be  emphasized  that  mjj  and  m^AT  are  normal  to  the  surface.  In  prac¬ 
tice,  these  quantities  are  likely  to  be  expressed  as  response  functions  defined 
for  the  limiting  cases.  For  purposes  of  illustration,  consider  pressure  coupling 
only,  end  define  R^fl,  Rbi^y 


=  Rbi  i 


all  parallel  acoustic  velocity  (A- 16) 


IK. 


_  _  A 

-3—  +  y  1  =  R.  no  parallel  acoustic  velocity  (A- 17) 

m,  a  X  p 


The  preceding  argument  shows  that  the  coupling  in  general  is  represented  by 

A  A 

Rfell  sin0  cosB  ,  (A-18) 

P  P 

and  the  rate  of  work  done  by  processes  at  the  element  of  area  dS  is 


[R^  sin20  +  R^^  cos20]  dS 


(A- 19) 


The  numerical  values  for  Rfcib  R^x  are  chose  either  calculated  or  measured  for 
the  appropriate  limiting  cases.  For  example,  data  taken  with  an  end-burning 
T-buraer  should  in  principle  provide  Rbi;  and  from  measurements  with  a  T-burner 
having  lateral  grains,  one  may  be  able  to  determine  R^j).  See  Section  9  for 
additional  discussions. 
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